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Abstract: The partition function of AA = 6 supersymmetric Chern-Simons-matter theory 

' ^ ' (known as ABJM theory) on S^, as well as certain Wilson loop observables, are captured by 

•'^ a zero dimensional super-matrix model. This super-matrix model is closely related to a matrix 

^L( model describing topological Chern-Simons theory on a lens space. We explore further these 

^ recent observations and extract more exact results in ABJM theory from the matrix model. In 

I— I particular we calculate the planar free energy, which matches at strong coupling the classical IIA 

-^ supergravity action on AdS4 x CP^ and gives the correct A'^" scaling for the number of degrees 

J> of freedom of the M2 brane theory. Furthermore we find contributions coming from world-sheet 

r instanton corrections in CP^. We also calculate non-planar corrections, both to the free energy 

QQ and to the Wilson loop expectation values. This matrix model appears also in the study of 

C^ topological strings on a toric Calabi-Yau manifold, and an intriguing connection arises between 

f — . the space of couplings of the planar ABJM theory and the moduli space of this Calabi-Yau. 

^^ In particular it suggests that, in addition to the usual perturbative and strong coupling (AdS) 

^^ expansions, a third natural expansion locus is the line where one of the two 't Hooft couplings 

• • vanishes and the other is finite. This is the conifold locus of the Calabi-Yau, and leads to an 

• ^H expansion around topological Chern-Simons theory. We present some explicit results for the 

^ partition function and Wilson loop observables around this locus. 
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1. Introduction and summary 

The discovery of Aharony, Bergman, Jafferis and Maldacena (ABJM) of the world- volume theory 
of coincident M2-branes [1] (following Bagger-Lambert and Gustavsson [2, 3]) provides a new 
interacting field theory with well defined weak and strong coupling expansions. A great deal 
of effort has been given to studying these two limits of the theory: three dimensional M = 6 
supersymmetric Chern-Simons-matter and type IIA string theory on AdS4 x CP^ (or M-theory 
on AdS4 X S'^/Zfc). For better or worse, both descriptions of the theory are much harder than the 
D3-brane analog: 4d A^ = 4 SYM and type IIB string theory on AdSs x S^. At weak coupling 
perturbative calculations in ABJM theory are rather subtle and for many quantities are in even 
powers of the coupling, while at strong coupling the geometry of CP'^ is more complicated than 
S^ and has, for example, non-trivial 2-cycles. 

An important breakthrough, which is the underpinning of the present study, was the work 
of Kapustin, Willett and Yaakov [4], who use the localization techniques of [5] to reduce the 
calculation of certain quantities in the gauge theory on S^ to finite dimensional matrix integrals^. 
These matrix integrals can be evaluated in a systematic expansion in 1/A^. Indeed, they have a 
natural supergroup structure, i.e., they are super-matrix models [7, 8], and are related to some 
previously studied bosonic matrix models [9, 10] by analytical continuation [8]. 

The solution of this matrix model allowed for the evaluation of the first exact interpolating 
function in this theory [8] giving a closed form expression for the expectation value of the 1/2 
BPS Wilson loop operator of [7] at all values of the coupling. This expression derived from 
the matrix-model reduction of the gauge theory reproduces exactly the known leading strong 
coupling result, the classical action of a macroscopic string in AdS4. 

The purpose of this paper is to explore further what can be learnt from the matrix model 
and its solution to the understanding of the physical 3d gauge theory and its string/M-theory 
dual. 

This is a broad subject, connected through the matrix model to special geometry, Chern- 
Simons (CS) theory, topological strings and more. One of the avenues we explore is the relation 
between the moduli space of the matrix model and the space of couplings of the gauge theory. It 
is very useful to consider the generalization of the gauge theory where the rank of the two gauge 
groups are not equal [11].^ The space of couplings is two dimensional and upon complexification, 
it matches the moduli space of the Riemann surface solving the planar matrix model. This 
surface is also the mirror to a well studied toric Calabi-Yau manifold known as local Fq, where 
Fq = P^ X P^ is a Hirzebruch surface. As we review in Section 3, this moduli space has three 
special loci: the orbifold point, the large radius limit and the conifold locus. 

These can be identified in the gauge theory respectively as the weakly coupled gauge theory, 
the strongly coupled theory described by string theory on AdS, and lastly the conifold locus is 
where the rank of one of the gauge group vanishes, so ABJM theory reduces to topological CS 
theory [12]. The first two are known duality frames with the AdS/CFT rules on how to evaluate 
observables on both sides. The simplicity of the conifold locus suggests that there should be 
another duality frame where ABJM theory is considered as a deformation of topological CS 
theory. We explore this in Section 6, where we calculate the partition function and Wilson 
loop observables around this point. It would be very interesting to learn how to calculate other 
quantities in this regime. 

^Similar results apply also to other 3d theories with J\f = 2 supersymmetry [4, 6]. 

^Though commonly known as ABJ theory, for simplicity we still call the theory with this extra parameter as 
ABJM theory. When specializing to the case of equal rank we refer to it as the "ABJM slice" . 



We present the matrix model for the ABJM theory and that for CS theory on the lens 
space L(2, 1) = S^/Z2 in the next section. The matrix model of ABJM has an underlying 
U{Ni\N2) symmetry while that of the lens space has U{Ni + A''2) symmetry, which in both cases 
are broken to U{Ni) x U{N2)- It is easy to see that the expressions for them are related by 
analytical continuation of A'^2 — ^ — ^^^2) or analogously a continuation of the 't Hooft coupling 
N2/k — )■ —N^jk (which may be attributed to the negative level of the CS coupling of this group 
in the ABJM theory). We can then go on to study the lens space model and analytically continue 
to ABJM at the end. 

Conveniently, the lens space matrix model has been studied in the past [10, 13, 14, 8]. The 
planar resolvent is known in closed form and the expressions for its periods are given as power 
series at special points in moduli space. We review the details of this matrix model and its 
solution in Sections 2 and 3. 

The matrix model of ABJM theory was derived by localization: it captures in a finite 
dimensional integral all observables of the full theory which preserve certain supercharges. At 
the time it was derived in [4], the only such observable (apart for the vacuum) was the 1/6 BPS 
Wilson loop constructed in [15, 16, 17] and 1/2 BPS vortex loop operators [18]. Indeed, the 
expectation value of the 1/6 BPS Wilson loop can be expressed as an observable in the ABJM 
matrix model, and by analytical continuation in the lens space model. 

Another class of Wilson loop operators, which preserve 1/2 of the supercharges, was con- 
structed in [7] and studied further in [19]. It is the dual of the most symmetric classical string 
solution in AdS4 x CP^. This Wilson loop is based on a super-connection in space-time and 
reduces upon localization to the trace of a supermatrix in the ABJM matrix model [7]. The 
different 1/2 BPS Wilson loops are classified by arbitrary representations of the supergroup 
L'^(A^i|A'^2)) and the 1/6 BPS ones are classified by a pair of representations^ of U{Ni) and U{N2) 
We will mostly concern ourselves with the 1/2 BPS Wilson loop in the fundamental represen- 
tation of U{Ni\N2) and the 1/6 BPS Wilson loop in the fundamental representation of U{Ni). 
The exception is Section 8.2 and Appendix B, where we study the 1/2 BPS Wilson loop in large 
symmetric and antisymmetric representations. There we also make contact with the vortex loop 
operators of [18]. 

Of course, the natural observables in CS theory are the partition function and Wilson loops, 
so these quantities were also studied earlier in the matrix models of CS (see, for example, [9, 10, 
20, 13, 21, 22]). This information is encoded in different period integrals on the surface solving 
the matrix model, as we explain in Section 2.2. It turns out that the 1/6 BPS loop is captured 
by a period integral around one of the two cuts in the planar solution and the 1/2 BPS Wilson 
loop by a period integral around both cuts, or alternatively, around the point at infinity, and is 
much easier to calculate [8]. 

With all this machinery presented in Sections 2 and 3 in hand, we are ready to calculate, 
and in Sections 4, 5 and 6 we study the partition function and Wilson loop observables in the 
three natural limits of the matrix model. First, in Sections 4 we look at the orbifold point, which 
is the weak coupling point of the matrix model and likewise of the physical ABJM theory. The 
calculations there are straight-forward and we present the answers to these quantities. A single 
term (1/6 BPS loop at 2-loops) was calculated independently directly in the field theory. All 
other terms are predictions for the higher order perturbative corrections. 

^Special combinations of representations of U{Ni) x U{N2) arc also representations of !7(A''i|A'^2), and in this 
case the 1/6 BPS and 1/2 BPS loops will have the same expression in the matrix model and the same VEVs. The 
proof of localization for the 1/2 BPS loop [7] relied on this equivalence. 



Section 5 addresses the strong coupling limit of the theory, where the matrix model should 
reproduce the semiclassical expansion of these observables in type IIA string theory on AdS4 x 
CP^. The expectation value of the Wilson loop was already derived in [8] and matched with a 
classical string in AdS. We first generalize the strong coupling expansion for the case of A^^i 7^ N2-, 
which corresponds to turning on a B-field in the AdS dual. This version of the theory was studied 
in [11] and a more precise analysis of the dictionary, capturing shifts in the charges, was presented 
in [23, 24]. Interestingly, it turns out that the matrix model knows about these shifted charges, 
and the strong coupling parameter turns out to be exactly the one calculated in [24] , rather than 
the naive coupling. 

In the same section we present also the calculation of the free energy in the matrix model. 
The result is proportional to iV^/yA (or a slight generalization for A'^i 7^ A''2). This scales at large 
N like A^^ " , which is indeed the M-theory prediction for the number of degrees of freedom on N 
coincident M2-branes [25]. Comparing with a supergravity calculation, we find precise agreement 
with the classical action of AdS4 x CP^. This is the first derivation of this large N scaling in the 
field theory side. The matrix model also provides an infinite series of instanton/anti-instanton 
corrections to both the partition function and to the Wilson loop expectation value, which we 
interpret as fundamental strings wrapping the CP^ inside CP^. 

We then turn to a third limit of the theory, when one of the gauge couplings is perturbative 
and the other one not. In the strict limit the ABJM theory reduces to topological CS and in the 
matrix model one cut is removed. We show how to perform explicit calculations in this regime 
both from the planar solution of the matrix model and directly by performing matrix integrals. 
In both approaches one can see the full lens space matrix model arising as a (rather complicated) 
observable in topological CS theory on S^. We speculate on possible tools of calculating directly 
in ABJM theory in this limit, where integrating out the bi- fundamental matter fields leads to 
correlation functions of Wilson loops in CS theory. We demonstrate the idea in the case of the 
1/6 BPS Wilson loop, which has a relatively simple perturbative expansion. This limit of the 
spin-chain of ABJM theory was considered in [26], and a similar system in four dimensions was 
studied in [27]. 

The brave souls that will make it to Sections 7 and 8 will find some new results on the 
non-planar corrections to the matrix model, and hence to ABJM theory. In section 7 we show 
that the full 1/A^ expansion of the free energy on S^ is completely determined by a recursive 
procedure based on direct integration [28, 29] of the holomorphic anomaly equations [30]. The 
ability to determine the full expansion is closely related to the integrability of topological string 
theory on toric Calabi-Yau threefolds (as discussed in for example [31]). By the AdS/CFT 
correspondence, the 1/A^ expansion obtained in this way determines the partition function of 
type IIA theory on the AdS4 x CP'^ background at all genera. This result is reminiscent of the 
"old" matrix models for non-critical strings, where a double-scaled 1/A^ expansion, encoded in 
an integrable system, captures the all-genus partition function of a string theory. The recursive 
procedure for the computation of the 1/A^ expansion is quite efficient in practice, and one can 
perform explicit computations at high genus. This allows us to study the large genus behavior 
of the 1/A^ corrections, and we check that they display the factorial growth ~ (25)! typical of 
string perturbation theory [32]. A careful examination of the coefficients suggests that this 1/A^ 
expansion is Borel summable. 

In Section 8 we present the genus one correction to the Wilson loop and expand it at both 
weak and strong coupling. Another topic covered there is that of "giant Wilson loops" [33, 34, 35], 
where in the supergravity dual (at least in AdSs x S^) a fundamental string is replaced by 
a D-brane. This happens for Wilson loops in representations of dimension comparable to N . 



We calculate the corresponding object in the matrix model and compare it to the vortex loop 
operators of [18]. 

One point we have not touched upon is the connection to topological strings. Since CS and 
the matrix model are related to topological strings, we expect there to be a direct connection 
between ABJM theory and a topological string theory. All the quantities captured by the matrix 
model should exist also in a topologically twisted version of ABJM theory, possibly along the 
lines of [36]. 

2. The ABJM matrix model and Wilson loops 

2.1 The matrix model and its planar limit 

The ABJM matrix model, obtained in [4], gives an explicit integral expression for the partition 
function of the ABJM theory on S^, as well as for Wilson loop VEVs. This matrix model is 
defined by the partition function 

Zabjm{Ni,N2, gs) 

. .}±^ ^fV d,. f^ d., n.<, (2sinh {^)y (2sinh (^))\ .^,^^.,,.^^.,.. 

.y 2-/^2. n.,(2coshr--.'^ 



(2.1) 
where the coupling gg is related to the Chern-Simons coupling k of the ABJM theory as 

2vri , , 

9s = -y. (2.2) 

In writing this matrix integral we have been very careful with its precise overall normalization, 
since one of our goals in the present paper is to compute the free energy on the sphere at strong 
coupling. The calculation of [4] captures the full k dependence of the partition function, but we 
have to fix an overall /c-independent normalization. This is done in two steps. First, we require 
that the above matrix integral reduces to the partition function for Chern-Simons theory on 
§^ when A'^i = or A^2 = (in a specific framing of S^). Once this is done, there is still a k- 
independent normalization factor which appears as a constant coefficient multiplying the cosh in 
the denominator. This term was not fixed in [4] , but it can be easily obtained from the formulae 
they presented. This calculation can be found in Appendix A, and leads to the matrix integral 
(2.1). 

The ABJM matrix model is closely related to the L(2, 1) lens space matrix model introduced 
in [9, 10]. This matrix model is defined by the partition function 
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(2.3) 
The relation between the partition functions is simply [8] 

ZABm{Ni,N2,gs) = Z^2,i){Ni,-N2,gs). (2.4) 



Since the large N expansion of the free energy gives a sequence of analytic functions of A'^i , N2 , 
once these functions are known in one model, they can be obtained in the other by the trivial 
change of sign A'^2 — ^ — -^''2. 
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Figure 1: Cuts in the z-plane and in the Z-plane. 

Let us now discuss the large N solution of the lens space matrix model, following [10, 13, 8]. 
At large A'', the two sets of eigenvalues, ^Uj, i^j, condense around two cuts. The cut of the in 
eigenvalues is centered around z = 0, while that of the Vi eigenvalues is centered around z = vri. 
We will write the cuts as 



in terms of the endpoints A, B. It is also useful to use the exponentiated variable 

Z = e^ 
In the Z plane the cuts (2.5) get mapped to 

(l/a,a), (-l/6,-&), 



a = e 



6 = e^, 



(2.5) 



(2.6) 



(2.7) 



which are centered around Z = 1, Z = — 1, respectively, see Fig. 1. We will use the same notation 
Ci^2 for the cuts in the Z plane. The large iV solution is encoded in the total resolvent of the 
matrix model, uiz). It is defined as [13] 

-(^)=5.(Tr(|^))=5.^E-th(^)^+5.^£tanh(^)^ (2-8) 



where 



U 



-^"1 



(2.9) 



We will denote by a;o(2) the planar limit of the resolvent, which was found in explicit form in 
[13]. It reads. 



ixiQ^z) = 2 log 



,-t/2 



where 



^/{Z + h){Z + 1/6) - ^{Z -a){Z- l/a) 



t = ti + t2 



(2.10) 
(2.11) 



is the total 't Hooft parameter. It is useful to introduce the variables 

/3 is related to the total 't Hooft parameter through 

/3 = e*. (2.13) 

All the relevant planar quantities can be expressed in terms of period integrals of the one-form 
LUo{z)dz. The 't Hooft parameters are given by 

ti = 2^i ojo{z)dz, i = l,2. (2.14) 

J Ci 
The planar free energy Fq satisfies the equation 

where the V cycle encloses, in the Z plane, the interval between —1/6 and 1/a (see Fig. 1).^ 

The derivatives of these periods can be calculated in closed form by adapting a trick from 
[37]. One finds. 



^'''' ' ^ ^^ ' ^"^ K{k), (2.16) 



dC 47ri / y/(Z2 - C^ + 1)2 - 4/^2^2 7r(l + ab) 

Cl,2 

and similarly 

1^ = -2^^?Hx (K{k) - :^ n(ni|fe) - -^ U{n,\k) ] , (2.17) 

op 7r(l + ab) \ I + ab I + ab 

where 

a + bV l-a2 6(a2-l) 



' 1 + abJ I + ab a[l + ab) 



Likewise for the period integral in (2.15) we find 



dC l + a6 ^ ^' 



ap 1 + ao \ [l + ab)[a + b) 



(2.19) 



where 

l + a6' ^ 6(l + a6) ^ 1 + a6 ^ ^ 

*The variable /? is related to the variable ^ in [8] hy (3 — ^/2. 

^Likewise one can calculate the second "B-cycle" period, and it will arise when solving the Picard-Fuchs 
equations at strong coupling in Section 3.2. 



We can now use the dictionary between the lens space matrix model and the ABJM matrix 
model given by (2.2) and (2.4) to get the planar solution of the latter model. In particular, the 
natural 't Hooft parameters in the ABJM theory 

A, = ^ (2.21) 

are obtained from the planar solution of the lens space matrix model by the replacement 

fi = 27riAi, t2 = -27riA2. (2.22) 

Since in the ABJM theory the couplings Ai,2 are real, the matrix model couplings ti,2 are pure 
imaginary. Thanks to (2.13) we know that /3 is of the form 

/3 = e2^i(^i-^2) (2.23) 

i.e., it must be a phase. 

For later convenience we introduce yet another parameterization of the couplings in terms 
of B and k 

B = Xi-X2 + ^, K = e--'^C- (2.24) 

B is identified as the B-field in the dual type IIA background [24]. Notice that it has a shift 
by —1/2 as compared to the original prescription in [11]. Clearly, all calculations in the matrix 
model are periodic under B — t- B + 1, up to possible monodromies (see (5.14) below). As we 
shall see later, the parameter k is real for physical values of Ai^2- 

2.2 Wilson loops 

One of the main results of [4] is that the VEV of the 1/6 BPS Wilson loop in ABJM theory, 
labelled by a representation R or U{Ni), can be obtained by calculating the VEV of the matrix 
e'^' in the matrix model (2.1), i.e., 

{W'J'') = gs (Tr^ (e'^0>ABJM mm , (2-25) 

1/2 

A 1/2 BPS loop W^ was constructed in [7] , where 7^ is a representation of the supergroup 
U{Ni\N2). in [7] it was also shown that it localizes to the matrix model correlator in the ABJM 
matrix model 

(W^i/')=5s(Str7eC/)ABJMMM> (2-26) 

with the same U as in (2.9). Though at first sight the minus sign on the lower block of U , may 
look surprising, it can be attributed to the fact that the Vj eigenvalues are shifted by vri from the 
real line. Due to the relation between the ABJM matrix model and the lens space matrix model, 
these correlators can be computed in the lens space matrix model as follows: 



(<^)=ff,(Tr«(e^0)L(2,i) 



^^^-^^ (2.27) 

l/2v 



N2->~N2 

where the super-representation TZ is regarded as a representation of U{Ni + A'^2)- 



To evaluate the Wilson loop one uses the resolvent, or equivalently, the eigenvalue densities 

(2.28) 



p«(Z)dZ = -^^ {^{Z + ie) - uj{Z - i6)) , Z G Ci, 



p(2)(z)dZ=^^(^(Z + i6)-a;(Z-ie)), Z ^ C^. 

which are each normalized in the planar approximation to unity 

/ pf AZ = 1. (2.29) 

For the 1/6 BPS Wilson loop in the fundamental representation one needs to integrate e^ = Z 
over the first cut 

Wll^^ = ti I p^^\Z)ZdZ= I ^uj{Z). (2.30) 

The correlator relevant for the 1/2 BPS Wilson loop (again in the fundamental representation) 
is much easier, since 

Wll^) = ti I p^^\z)ZdZ-t2 [ p^^'>{Z)ZdZ= I ^uj{Z) (2.31) 

/ Jci Jc2 Joo 47ri 

and it can be obtained by expanding uj{Z) around Z — )■ oo. 

The comparison to the case of the 1/2 BPS Wilson loop in AA = 4 SYM in 4d is straight- 
forward. In that case the matrix model is Gaussian and the eigenvalue density in the planar 
approximation follows Wigner's semi-circle law. Doing the integral with the insertion of e^ gives 
a modified Bessel function [38] 

P0{Z) = -^VX-Z^ => (W^id'Ar = 4)planar= / Po{z) e' dz = ^h{^). (2.32) 

For the ABJM matrix model all the expressions are more complicated. Still the derivative 
with respect to ( and /? of the integral expression for the 1/6 BPS Wilson loop (2.30) can be 
written in closed form [8], like the integrals (2.16) and (2.17) 

d^{W^^^) = -^ ^] {aK{k) - (a + b)U{n\k)) 



TT Vab(l + ab) 

^ ' (2.33) 



For the 1/2 BPS Wilson loop of [7] the situation is much simpler and in the planar approximation 
one needs only the large Z behavior of uq (2.10) 



One finds [8] 



^ -, I C I C^ + 2/3^-2 C(C^ + 6/3^-3) 4 

^o-t+- + —, + —, + 0{Z ). 



.l/2v _ C 



(2.34) 



(^^□'Oplanar = \, (2.35) 



which can then be expanded in different regimes. We will elaborate on the expansion of this 
expression in the next sections and will also turn to the non-planar corrections to it and to that 
of the 1/6 BPS loop in Section 8. 

As a simple generalization, by the replacement Z — )■ Z' on the right hand side of (2.31), the 
higher order terms in the expansion (2.34) give the expectation values of multiply wrapped 1/2 
BPS Wilson loops where [/ — t- [/' in (2.26). For even winding the sign in the lower block of the 
matrix U (2.9) is absent. This is consistent with the gauge theory calculation [7], where this 
sign arose from the requirement of supersymmetry invariance in the presence of the fermionic 
couplings which are antiperiodic, as should be the case for a singly-wound contractible cycle (see 
also the discussion in [19]). 

The normalization of the Wilson loop as given by (2.30) and (2.31) is not the same as in the 
4d AA = 4 case (2.32). For the 1/6 BPS loop, the leading term at weak coupling is ti = 2TriNi/k. 
This means that the trace in the fundamental is normalized by a factor of gs- For the 1/2 BPS 
loop the leading term is ti zt t2 = dsi^^i T -^2), where the sign depends on the winding number. 
We will comment more about this normalization in Section 5.2. 

3. Moduli space, Picard Fuchs equations and periods 

In this section we present the tools for solving the lens space matrix model using special geometry. 
We present three special points in the moduli space of the theory and write explicit expressions 
for the four periods of ojq at the vicinity of these points. 

The lens space matrix model is equivalent to topological string theory on local Fq = P""^ x P^. 
The 1/A^ expansions of the free energy and of the 1/2 BPS Wilson loop VEV are the genus 
expansions of closed and open topological string amplitudes. The planar content of the theory is 
encoded in the periods of the mirror geometry described by the family of elliptic curves S, which 
can be written as 

zix^ + a; + 1 - a/(1 + X + z\x^)'^ - ^z^x^ ,^^, 

y = 2 • ^ > 

Here, z\^zi parametrize the moduli space of complex structures, which is the mirror to the 
enlarged Kahler moduli space of local Fq. This moduli space has a very rich structure first 
uncovered in [10] and further studied in, for example, [37, 31] by using the standard techniques 
of mirror symmetry. 

Notice that the mirror geometry (3.1) is closely related to the resolvent u;q{Z). Indeed, one 
finds that ujq{Z) '^ log y{x) provided we identify the variables as 

(3.2) 
(3.3) 



X = - 


-Zz^ 




and 
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This can also be expressed as (2.24) 
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(3.4) 

Let us now discuss in some detail the moduli space of (3.1), since it will play a fundamental 
role in the following. It has complex dimension two, corresponding to the two complexified Kahler 
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parameters of local Fq. The coordinates zi, Z2 (or ^, /3) are global coordinates in this moduli space. 
Another way of parametrizing it is to use the periods of the meromorphic one-form 

uj = logy{x) — (3.5) 

X 

As it is well-known, these periods are annihilated by a pair of differential operators called Picard- 
Fuchs operators. In terms of zi,Z2, the operators are 

£l = Z2{1 - 4z2)ei - 4z?ef - 8ziZ2^iC2 - Qzi^i + (1 - 6Z2)6, 

(3.6) 
£2 = zi{l - Azi)f^ - Azlil - 8ziZ2Cia - QZ2C2 + (1 - 6zi)6, 

where 

«• = !-■ <^^^' 

These operators lead to a system of differential equations known as Picard-Fuchs (PF) equations. 
An important property of the moduli space is the existence of special points, generalizing the 
regular singular points of ODEs on C. The PF system can be solved around these points, and 
the solutions give a basis for the periods of the meromorphic one-form. We can use two of 
the solutions to parametrize the moduli space near a singular point, and the resulting local 
coordinates, given by periods, are usually called flat coordinates. 

3.1 Orbifold point, or weak coupling 

There are three types of special points in the moduli space. The first one is the orbifold point 
discovered in [10], which is the relevant one in order to make contact with the matrix model. To 
study this point one has to use the global variables 

xi = 1 , X2 = 7 ^. (3.8) 

The orbifold point is then defined as xi = X2 = 0, and in terms of these variables the Picard- 
Fuchs system is given by the two operators 

Ci=-{8- 8x1 + xl)x2dx2 - - (4 - (2 - xifxl) dl^ - xi(2 - 3xi + x\)x2dx^dx2 
- (1 - xi)xldxi + (1 - xifxldl^, 

£2 = (2 - xi)x25:r2 - (1 - (1 - xi)xl)dl^ - x\dx^ - 2(1 - xi)xiX2dxidx2 + (1 - xi)xldl^. 

(3.9) 
A basis of periods near the orbifold point was found in [10]. It reads, 

CTi = -log(l -xi), 

^2 = / Cm,nX\ X2 , 

^ (3.10) 



-7^2 = 0"2 log Xl + ^ dm,nX'^X2 
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where the coefficients Cm,n and dm,n vanish for non-positive n or m as well as for all even n. They 
satisfy the following recursion relations with the seed values ci^i = 1, di^i = and di^s = —1/6: 



Cm,n — ^/-___ __\/-___ -,\ Cm— l,n) 



Cm,n— 2) 



n + 2 - 2m)2 
A{m — n){m — 1) 

(n — 2)^(m — n + 2)(m — n + 1) 
n(n — l)(2r?T, — n)2 

_ (n + 2 - 2m)^dm-i,n + 4(n^ - n - 2m + 2)cm^„ 
4(?TT, — 1)(?TT- — n)(n + 2 — 2?7i) 

_ (n-2)2(m-n + l)(m-n + 2) /I 1 4 \ 

dm,n — 7 TTTo V? "m,n-2 + I "T^ H ~~r H ^^ I C-mn- 

n{n — \)[lra — n)'^ \m — n + Z m — n + 1 n — Ira J 

(3.11) 
The 't Hooft parameters of the matrix model are period integrals of the meromorphic one-form, 
therefore they must be linear combinations of the periods above, and one finds [10] 

h = ^{(Ti+a2), t2 = -{(Ti-a2). (3.12) 

An expansion around the orbifold point leads to a regime in which ti,t2 are very small. In view 
of (2.22) this corresponds, in the ABJM model, to the weakly coupled theory 

Ai, A2 < 1. (3.13) 

The remaining period in (3.10) might be used to compute the genus zero free energy of the matrix 
model. Using the normalization appropriate for the ABJM matrix model, we find 

3.2 Large radius, or strong coupling 

The second point that we will be interested in is the so-called large radius point corresponding 
to zi = Z2 = 0. This is the point where the Calabi-Yau manifold is in its geometric phase, and 
the expansion of the genus zero free energy near that point leads to the counting of holomorphic 
curves with Gromov-Witten invariants. The solutions to the Picard-Fuchs equations (3.6) near 
this point can be obtained in a systematic way by considering the so-called fundamental period 

wq{zi,Z2;pi,P2)= > -^^^ — ^o \r^^^u^ — \2r(^ ^/^ — yT'^i ^2 • 3.15 

kJ>Q ^^1 + 2/92)r(l + k + piY T{l + l + piY 

As reviewed in for example [39], a basis of solutions to the PF equations can be obtained by 
acting on the fundamental period with the following differential operators 

Df^=d,^, Df^ = \K,,udp,dp,. (3.16) 

Here Kijk are the classical triple intersection numbers of the Calabi-Yau. This leads to the 
periods 

Ti{zi,Z2) = -Dl 'zuo{zi,Z2;pi,P2) 

''=''=' (3.17) 

(2) 

Fi{zi,Z2) = Dl 'wo{zi,Z2;pi,P2) 

Pl=P2=0 
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These periods should be hnearly related to those defined in the matrix model in equations (2.14) 
and (2.15). We present now some explicit expressions for them that we will use in Sections 5.1 
and 5.3 to solve for these relations (see equations (5.3) and (5.23)). 

In general, one normalizes these periods and divides them by the fundamental period eval- 
uated at pi = p2 = 0. But in local mirror symmetry we have [40] 



-^Qizi, Z2; pi, P2] 



1. 



P1=P2=0 



(3.18) 



The Ti are single-logarithm solutions, and they are identified in standard mirror symmetry with 
the complexified Kahler parameters, while the Fi are double- logarithm solutions and they are 
identified with the derivatives of the large radius genus zero free energy w.r.t. the Tj. In our 
case, we find the explicit expressions 



-Ti =logzi +a;(^) (2:1,2:2), 
-T2 =log^2 +0;^^) (2:1, 2:2), 



(3.19) 



where 



a;W(zi,Z2) = 2 Y. m 



T{2k + 2l) 



k,l>0, 
(k,l)^(0,0) 



(i + fc)2r(i + /)2 



^5^4 = 2zi + 2z2 + 34 + 12^122 + 3z^ + --- (3.20) 



In order to obtain the Fi we have to compute the double derivatives w.r.t. the parameters pi, 
P2- We find 



5pi^o(^i,^2;/>i,P2) 



,(1) 



(2), 



P1=P2=0 



log Zi + 2logZiUJ^ '{Z1,Z2) +UJi {zi,Z2) 



where 



Similarly, 



where 

,(2), 



f\zuZ2) = ^ Y, ^^^f^w^^! n2 (^(2fc + 20 - ^(1 + k)) z\z\. 



k,l>0, 
(k,l)jt(0,0) 



+ A;)2r(l + /)2 



dp^'cuo{zi,Z2;pi,P2) 



Pl=P2=0 



log 2:2 + 2 log 2:2a;'' -'(21,^2) + ^2 izi,Z2 



m, 



(3.21) 
(3.22) 

(3.23) 



4^)(.i,.2) = 8 Y. fa 
Finally, 



T{2k + 2l) 



k.l>0, 
(fe, 07^(0,0) 



+ A:)2r(l + /)^ 



(^(2A: + 21) - V(l + /)) ^i4 = 0^1 {Z2, zi). (3.24) 



,(2)/ 



dp,dp2'!Uo{zi,Z2;pi,P2) 



pl=P2=0 



log 2:1 log 22 + (log 2:1 +log22)w^ 'izi,Z2 
+ 2 [^l (zi, Z2) + U;2 {Zi, Z2)j . 



(3.25) 



The double log periods are obtained as linear combinations of the above, by using the explicit 
expressions for the classical intersection numbers that can be found in for example [31] 



1 



1 



Kill 



K112 



K122 



1 

'4' 



1 



K222 



(3.26) 
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We find: 



' "i2 , ,_ o ri , ._ r-i2 



^^1(21,2:2) = -- (^pi^o - 2Dp^p2UJo - -Dpi wo) 

= -- (log^Zl - 21og2;ilogZ2 - log^ Z2) + -logZ2Uj''^\zi,Z2) + -UJ2 (^1)^2), 
o 4 o 

F2{zi, Z2) = -- (-Dpj^o - 2Dp^p^u}Q - DI^ujq) 

= -- (-log^Zl - 21ogZllogZ2 + log^Z2) + -\0gZlUJ^'^\zi,Z2) + -uf\zi,Z2). 
o 4 o 

(3.27) 

They satisfy the symmetry property 

Fi{zi,Z2) = F2{z2,zi). (3.28) 

The reason why we are interested in the large radius point is because it describes the structure 
of the AB JM theory at strong coupling. In the region where Z2 is small, X2 is large and the periods 
^1,2 grow. In general, the expansions of the periods around the special points have a finite radius 
of convergence, but they can be analytically continued to the other "patches" . Since their analytic 
continuation satisfies the PF equation, we know for example that the analytic continuation of 
the orbifold periods to the large radius patch must be linear combinations of the periods at large 
radius. This provides an easy way to perform the analytic continuation which will be carried out 
in detail in the Section 5, where we will verify that indeed the region near the large radius point 
corresponds to 

Ai, A2 > 1. (3.29) 

3.3 Conifold locus 

Finally, the third set of special points is the conifold locus. This is defined by A = 0, where 

A = 1 - 8(zi + Z2) + 16(zi - Z2f. (3.30) 

In terms of the variables C,, (3, this locus corresponds to the four lines 

C = -2/3±2, C = 2/3±2. (3.31) 

The conifold locus is the place where cycles in the geometry collapse to zero size. The first 
two lines correspond to a = ±1, i.e., the collapse of the Ci cycle, while the second set of lines 
corresponds to 6 = =Fl) ^-6., to the collapse of the C2 cycle. In principle we can solve the PF 
system near any point in the conifold locus, but in practice it is useful to focus on the point 

zi = Z2 = — (3.32) 

lb 

which has been studied in [31]. We will call it the symmetric conifold point. Appropriate global 
coordinates around this point are^ 

yi = 1 - -, y2 = 1 - T^. (3.33) 

Z2 IQzi 



^These are slightly different from the ones used in [31]. 
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In terms of these coordinates, the PF system reads 



A =dy, - 2(1 - y2)dl^ - 8(1 - yifdy, + 8(1 - yifdl^, 

£2 = - (7 - ^y2)dy, + 2(3 - 77/2 + 4yi)92^ - 8(1 - yi)dy, (3.34) 

2fl2 



- 16(1 - yi)(l - y2)dy,dy, + 8(1 - yifdy^. 

Notice that, strictly speaking, the orbifold point does not belong to the conifold locus, once the 
moduli space is compactified and resolved [10]. A generic point in the conifold locus has then 
ti = or ^2 = 0, but not both, and expanding around the conifold locus means, in the ABJM 
theory, an expansion in the region 

Ai < 1, A2 ~ 1, (3.35) 

or in the region with A2 exchanged with Ai. This regime of the ABJM theory has been considered 
in [26]. 

It was observed in [41] that the moduli space of the local Fq surface can be mapped to a well- 
known moduli space, namely the Seiberg-Witten (SW) ?i-plane [42]. This plane is parametrized 
by a single complex variable u. The relation between the moduli is 

^=^(/5 + r')-^. (3.36) 

The three singular points that we have discussed (large radius, orbifold, and symmetric conifold) 
map to the points u = 00, +1, —1. These are the semiclassical, monopole and dyon points of SW 
theory. As we will see, they can be identified with interesting points in ABJM theory. 

An important set of quantities in the study of moduli spaces of CY threefolds are the three- 
point couplings or Yukawa couplings, CziZjZk- These are the components of a completely symmet- 
ric degree three covariant tensor on the moduli space. When expressed in terms of flat coordinates 
they give the third derivatives of the genus zero free energy. In terms of the coordinates zi,Z2, 
the Yukawa couplings are given by [10, 31] 



Ci 



11 



<^112 - J^ 

C122 



(1 - 4^2)^- 16-^1 (1 + zi) 
AzfA 

I6z1 - (1 - 4^2)^ 



^ (3.37) 

16z|-(l-4zi)2 



G 



222 



Aziz^A 

{I - Azi)"^ - I6z2{l + Z2) 
4zfA 



3.4 The moduli space of the ABJM theory 

The matrix model of ABJM is closely related to the lens space matrix model, and therefore so are 
also the moduli spaces of the theories. Some of the explicit relations needed for this identification 
will be presented only in the following sections, but we would still like to present here the main 
points on the moduli space. 

We can think about the moduli space of the planar ABJM theory as the space of admissible 
values of the 't Hooft parameters Ai, A2. We will assume for simplicity that A: > 0. The theory 
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with negative values of k can be obtained from this one by a parity transformation. In the gauge 
theory Ai^2 must be rational and non-negative (for k > 0). Moreover, according to [11], any value 
of Ai^2 is admissible as long as 

|Ai-A2|<l. (3.38) 

This moduli space can be parametrized by the B field and k, which from the explicit expressions 
derived below (4.1) and (5.11) has to be real and positive. It can be identified as a real submani- 
fold of the moduli space of local Fq. Moreover, we can identify the singular points of this moduli 
space with natural limits of ABJM theory (see Fig. 2): 



X2 + 



21 , 22 — > 00 




Xi, X2 



B^O 



2 



B^ 1 



orbifold 



A 



1 



Figure 2: The moduli space of the ABJM theory, describing the possible values of the 't Hooft couplings 
Ai.2, can be parametrized by a real submanifold of the moduli space of local Fg, here depicted as a sphere. 
The orbifold point maps to the origin, while the conifold locus (which is represented by a dashed line) 
maps to the two axes. 



1. The weak coupling regime Ai^2 — ^ corresponds to the orbifold point of the local Fq geometry 
K = 0, B = 1/2. In terms of type IIA theory, this is also an orbifold geometry with a small 



radius but a nonzero value for the B field. 



2. The strong coupling regime Ai^2 
limit of the local Fq geometry. 



oo (where also k — ?■ oo) corresponds to the large radius 



3. Out of the four lines (3.31) in the conifold locus A = 0, only two lead to k € M. They are 
the curves in the (k, B) plane with k = ib4cos7ri?, which correspond respectively to a = 1 
and 6=1, therefore to Ai = or A2 = 0. Hence, the boundary of the ABJM moduli space 
given by min(Ai, A2) = corresponds to 



k{B) 



—icosnB . 
4 cos ttB , 



B> 1/2 
B < 1/2 



(3.39) 



In particular, the symmetric conifold point zi = Z2 = 1/16 corresponds to i? = n G Z, 
K = ±4. Along the curve (3.39), one of the two gauge groups of the ABJM theory is absent, 
so the theory reduces to topological CS theory. We examine this regime in Section 6. 
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Given a fixed value of the B field, we can describe the real one-dimensional moduli space of 

the ABJM theory as a real submanifold of the n-plane of Seiberg-Witten theory, by using (3.36) 

in the form 

^2 



cos(27rS) + 



(3.40) 



Singular points in moduli space become then the well-known singularities of SW theory. For 
example, when B = 1/2, the moduli space, described by «; G [0, cxd), becomes the region u E 
[1, oo). The orbifold point (weakly coupled ABJM theory) maps to the monopole point, while the 
large radius point (strongly coupled ABJM theory) corresponds to the semi-classical region (see 
Fig. 3). Notice that the conifold point would map to the dyon point of Seiberg-Witten theory, 
but this does not belong to the moduli space of ABJM theory with B = 1/2. We can however 
realize it by making an analytic continuation of the 't Hooft coupling to complex values. The 
dyon point corresponds then to the point k^ = —16, which leads by (5.5) to an imaginary value 



A 



2iK 



n^ 



(3.41) 



where K is Catalan's number. 

As usual, string dualities lead to a full complexification of the moduli space of 't Hooft 
parameters. In the case of ABJM theory, the complexified moduli space for the variables Ai^2 is 
simply the moduli space of the parameters /3, C; which is a Z2 x Z2 covering of the moduli space 
parametrized by 2:1,2. 



u = -l 



u = l 

• — 



u 



00 



conifold point 



weakly coupled 
ABJM theory 



strongly coupled 
ABJM theory 



Figure 3: The moduli space of the ABJM theory for B = 1/2 can be mapped to the line [l,oo) in the 
u plane of Seiberg-Witten theory, which is here shown in red. The monopole point corresponds to the 
weakly coupled ABJM theory, while the semiclassical limit corresponds to the strongly coupled theory. 



4. Weak coupling 

In principle, to study the matrix model at weak coupling one does not need the sophisticated 
tools presented in the previous section. One can do perturbative calculations directly in the 
integral expressions (2.1) or (2.3) for the matrix model. A calculation of the 1/6 BPS Wilson 
loop to three loop order was indeed done in this way in the original paper [4]. 
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Still, the explicit expressions for the periods ai^2 (3.10) and their relation to ti^2 (3-12) gives 
a much more efficient way to obtain perturbative, planar expansions. Inverting these relations 
we find the weak coupling expression for k (2.24) 



2i(ii - *2) - Y^ {4 + 3t?t2 - 3titi - 4) 

^ (tf + 5tit2 - wtl4 + Wtl4 - 5h4 - 4) + Oif). 



(4.1) 



This agrees with the weak coupling expansion of the inverse of the exact mirror map (5.5), 
obtained in [8]. 

Using the dictionary relating the 't Hooft couplings (2.22) we immediately get the result for 
the 1/2 BPS Wilson loop in the planar approximation (2.35) 



(Wll^) = e"'-^- =e"'(^i-^2) 2vri(Ai + A. 



° ' " 2 

^4 



2 



l-^(A?-4AiA2 + A^) 



,>4 
120 



(4.2) 
{Xf - 6A?A2 - 4A?Ai - 6AiAi + A^) + 0{X^ 



In this expression we factored out the term 27ri(Ai + A2), which depends on the overall nor- 
malization of the Wilson loop, as mentioned after (2.35). There is also the extra phase factor, 
which appears also at strong coupling and can be attributed to framing. Note that so far this 
expansion has not been reproduced directly in the gauge theory, as even the two-loop graphs are 
quite subtle. 

For the 1/6 BPS Wilson loop, using the explicit expression (2.33) and expanding at low 
orders one finds [8] 

^l/6^^ = e-^^27riAi (l - ^Ai(Ai - 6A2) - ^AiA^ + ^Ai (A? - 10A?A2 - 20A3) + 0{X'')] . 

(4.3) 
Again the exponent is a framing factor and the factor of 27riAi is due to the normalization chosen 
in (2.30). This expression agrees with the 2-loop calculations in [15, 16, 17]. Note that the 3-loop 
analysis in [17], done for Ai = A2, misses the next term, due to a projection which essentially 
removes all terms at odd orders in perturbation theory. 

Next we turn to the free energy. Here we notice that the period in (2.15) gives only the 
derivative of the free energy. Indeed, within the formalism of special geometry developed above, 
the planar free energy of the matrix model is only determined up to quadratic terms in the 
't Hooft couplings. These have to be fixed by direct calculation in the matrix model 

F = f log (^) + f log (^) - \(N^, + iVj) - l„g(4,Ar.». + . ^ ^ (4.4, 

The last term comes from the normalization of the cosh term in (2.1), while the remaining terms 
are just the free energies for two Gaussian matrix models with couplings ib27ri/fc. Notice that 
the above free energy has an imaginary piece given by 

'^{Ni-N2){{Ni-N2f-l). (4.5) 
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Using the identification of the periods at weak couphng (3.14) we write down the next term in 
the perturbative expansion 



72fc2 



{Nf - 6NfN2 + ISNfN^ - 6iVi iV| + iV|) . (4.6) 



It would be interesting to try to reproduce these expressions directly from studying perturbative 
ABJM theory on S^. 

5. Strong coupling expansion and the AdS dual 

We turn now to the strong coupling limit of the matrix model, where we have to find the analytic 
continuation of the 't Hooft parameters to the strong coupling region, as functions of the global 
parameters of moduli space. We will see how the shift of the charges discussed in [23, 24] emerges 
naturally from our computation. We will also evaluate the free energy in this regime and compare 
with the classical action of the vacuum AdS dual, deriving in this way the A^^' ^ behavior of the 
degrees of freedom. 

5.1 Analytic continuation and shifted charges 

In order to perform the analytic continuation of the 't Hooft parameters, we use the explicit 
representation of the periods in terms of integrals given in (2.14) as well as their derivatives 
(2.16)-(2.17). Let us start by discussing ti. We study its behavior at large ( but fixed /3, which 
is the large radius region. We find 

t = i'°<-f)+*"'' | = -2^('°«<-^^^'+-)+"a). (5.1) 

and this gives the leading behavior 

ti = -^(log(-C')+vri)log^ + ... (5.2) 

In the physical theory ti should be imaginary and /3 a phase. By examining (5.2), this implies 
that K is real. From (3.4) we then see that zi = Z2 and henceforth we label it zi = z. 

We know also that ti must be a linear combination of the periods at large radius. Using 
that zi = l/C,"^ and Z2 = {/3/Q'^, and comparing (5.2) to the behavior of the periods (3.19) and 
(3.27), we find 

h = —(Fi + F2) - -T2 - — , 

'-' ', " . (5.3) 

The constants ±111/6 cannot be fixed by using the above information, but they can be fixed by 
specializing to the ABJM slice zi = Z2, as we will see in a moment. 
A simple calculation leads to the following explicit expression 

^■<- «) 4 (^= - i) - ^ - ^ - 'W^'" (- ^^ - rb i^?' - -?') (- ^-' ■ <^^^' 

This expansion is valid in the region k — )■ +00. Notice that it is manifestly real when k is real 
and positive. 



19 



As a check of the above expression, we can particularize to the ABJM shce Ai = A2 = A, 
(B = 1/2), which corresponds in the gauge theory, to having identical gauge groups in the two 
nodes of the quiver, i.e., Ni = N2. The mirror map for this case was obtained in [8] as 

/ 1\ K /I 1 1 3 k2\ ^ ^ 

The strong coupling expansion of this expression at k S> 1 is 

^(-,B=l] = '^ + l- + Oi.~% (5.6) 



2j 27r2 24 

in agreement with (5.4). This also fixes the constants in (5.3). 

As in [8], the observables of the model are naturally functions of (", /3 (alternatively k, B), 
and we have to re-express them in terms of Ai^2- Equation (5.4) shows that the natural variable 
at strong coupling is not Ai, but rather 

A = A. - 1 (b^ - 1) - 1 = 1(A, + A,) - \(M - ^.f - ^. (5^T) 

In particular, it is only when expressed in terms of this variable that k is a periodic function of 
X,B. 

Remarkably, the above shift is precisely the one found in [24]. In the type IIA realization of 
the ABJ theory U{M2)k x U{M2 + M4)_fc, where M2 corresponds to the number of D2 branes 
and M4 to the number of D4 branes, the Maxwell charge of the D2 branes is not M2, but rather 

% = M.-^(B^4)-i(.-i). (5.8, 

where 

M4 1 , , 

After dividing by k and taking the large k limit, we recover (5.7) with 

A=^. (5.10) 

k 

The relation between A and k can be inverted at strong coupling, generalizing [8] to B ^ g) 
and it is of the form 

K(A,S) =e-^ (l + Y^ce l^,^] e-''-^] (5.11) 

\ ^ VvrV2A / / 

where 

2£-l 

q(x,/3) = 5^4')(/3)x^ (5.12) 

k=0 

The coefficients c^ (/3) are Laurent polynomials in j3,/3~^, of degree i, and symmetric under the 
exchange /3 -H- /?~^. In other words, they can be written as polynomials in cos(27rmi?), so they 
are periodic in B, with period 1. We find, for example, 



ci(x,/3) = -(/? + r')(i-|) 






■^ /^ --,0 ^ .^ ^ o\ OX / „ .-. 1 \ 2 X 
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The fact that q(x, /3) are polynomials in x of degree 2i — 1, rather than power series, comes out 
from an explicit calculation of the first few cases, and we have not established it. 

From the explicit expression (5.4) we can implement the symmetries of the model as a 
function of k and B (or equivalently, zi and ^2)- For example, the transformation 

Ni^2Ni + k- N2, N2 -^ Ni (5.14) 

simply corresponds to periodicity in the B field 

B^B + 1 (5.15) 

while K remains unchanged. From the point of view of the zi^2 variables, this is simply a mon- 
odromy transformation zi^2 — ^ e^^'^'zi^2- Notice that not all the values of k lead to admissible 
values of Ai^2) since min(Ai,A2) > 0. This means that the boundary of moduli space is the 
conifold locus (3.39). 

5.2 Wilson loops at strong coupling and semi classical strings 

As an application of the explicit expression for k (5.11), we can use (2.24) to immediately obtain 
the VEV of the 1/2 BPS Wilson loop (2.35) at strong couphng 

{W'j'),=0 = le-'^KCX,B). (5.16) 

Note that this is a real function of A, B, up to the overall phase involving the B field. This is 
the same phase that appears also in the weak-coupling result (4.2) and arises also in field theory 
calculations as a framing-dependant term [12, 43, 44]. The matrix model always gives the answer 
for framing=l. 

The result for the 1/6 BPS Wilson loop is, as usual more complicated, but can still be written 
in a power series expansion at strong coupling. We quote only the leading strong coupling result 
for Ai = A2 [8] 

{w'J'),=o^~e-''^e-^\ (5.17) 

We would like to comment about the normalization of the operators. As mentioned after 
(2.35), the normalization chosen there is such that the trace of the identity in the fundamental of 
U{Ni) gives ti = 2TT\Ni/k and for the fundamental of U{Ni\N2) (with a minus sign as in (2.9), 
it gives ti—t2 = 2it[{Ni + N2)/k. In CS theory these normalizations are quite common, but they 
may be not the most natural ones in the ABJM theory. 

An alternative normalization is to divide by this term, such that at weak coupling the 
expansion of the Wilson loop will be {W) ~ 1 + • • • . This is the normalization chosen in [8], and 
hence the slight differences in the preceding equations from that reference. Note, though, that 
with such a normalization, one would have to divide the doubly- wound 1/2 BPS Wilson loop in 
the fundamental representation by the super-trace of the identity, which is 27ri(A'^i — N2)/k and 
is singular for Ni = N2- 

There should be a natural choice of normalization that would reproduce the correct nor- 
malization fo the one- loop partition function of the classical string in AdS4 x CP^. To this day, 
though, a fully satisfactory calculation for the analog string in AdSs x S^ giving the factor of 
^-3/4 (jgrived from the the Gaussian matrix model does not exist. One argument, based on 
world-sheet arguments was given in [45], but it is not clear why this argument would be modified 
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for ABJM theory. Direct calculations of the determinant [46, 47] were not conclusive. A possible 
trick to derive it was proposed in [48] by considering a 1/4 BPS generalization of the circular 
Wilson loop, where three zero modes of the the Wilson loop of [49] are explicitly broken and the 
integral over them gives this factor. It would be interesting to construct such generalization to 
the Wilson loop of [7] and see if a similar argument can be derived from that. 

Regardless of the overall normalization, one can compare those of the 1/2 BPS loop and the 
1/4 BPS loop. Ignoring numerical constants and the framing factor, the ratio is 






A, (5.18) 



which is proportional to the volume of a CP^ inside CP^. Indeed, it was argued in [15, 17] that 
the string description of the 1/6 BPS Wilson loop should be in terms of a string smeared over 
such a cycle. 

5.3 The planar free energy and a derivation of the N'^''^ behaviour 

In this section we study the free energy at strong coupling. We derive the N^'"^ behavior charac- 
teristic of M2 branes [25] , and we match the exact coefficient with a gravity calculation in type 
II A superstring on AdS4 x CP^. 

The free energy of the matrix model has a large N expansion of the form 

oo 

F = logZ = J]<7'^-2f,(Ai, A2). (5.19) 

9=0 

This is the way the genus expansion is typically expressed in topological string theory. To 
compare with the gauge theory and the AdS dual one may choose to rewrite this series as an 
expansion in powers of 1/A^ by absorbing factors of A into Fg. 

As mentioned in Section 4, the formalism of special geometry determines the planar free 
energy only up to quadratic terms in the 't Hooft couplings, and these have to be fixed from the 
explicit weak coupling calculation in the matrix model (4.4). 

Let us now consider the derivative of the genus zero free energy (2.15), and study its analytic 
continuation to strong coupling as we have done for ti at the top of Section 5.1. Expanding (2.19) 
for large k we find 

f = -f + o(r^), | = ^(rt (5-20) 

so 

J= -vri log C + 0(C°) = -TTi log K + 7r2B + 0(^0, S°), k ^ 00, (5.21) 

From this leading large k behavior we have that in the ABJM slice 

2tt^V2X, (5.22) 



9Fo ^ ,„3 



dX 

which can be integrated to give the leading term in (5.34) and the match with the supergravity 
calculation presented below. 

But to get the full series of corrections we should proceed more carefully. We know that the 
result of the continuation should be a linear combination of periods, and comparing to (3.19) we 
see that we can express the period as 

_ TTit dFn dFt) 7ri ,„ „ , , , 
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The constant term can be fixed by looking at the solution on the ABJM slice A'^i = N2, which 
can be obtained as follows. Since on the slice we effectively have a one-parameter model, there 
is only one Yukawa coupling, which we can integrate to obtain Fq. From (3.37) we easily obtain 



dlFoiX) = ]cxxx 



1287r6 



1 



Ai=-A2 k{k'^ + 16) X (i'^Y 



(5.24) 



where the factor of 4 is introduced to match the normalization of the matrix model, and we used 

that 

dA 1 ^, /Ik' 



Integrating once, we find 



dK Att"^ \ 4 



K' ('—) 



where ai is an integration constant and we have used the Legendre relation 

E'K + EK' -KK' = -. 



(5.25) 
(5.26) 

(5.27) 



A further integration leads to the following expression in terms of a Meijer function 



%w) = fGS(y;4 



16 



-l-aiA-l- a2. 



(5.28) 



Comparison with the matrix model free energy at weak coupling (4.4) fixes oi = 47r^i, a2 = 0, 
so we can write 



%W) = fGg(y;j. 



K^\ TT^iK 



16/ "^ 2 



A^llAa' "'^ 



2'2'2' ' 2' 16 
If we integrate this expression with the following choice of integration constant. 



(5.29) 



Fo(A)= / dA'avFo(A') 
Jo 

we obtain the correct weak coupling expansion. 

We can now analytically continue the r.h.s. of (5.29) to k = 00, and we obtain 



(5.30) 



4tt 



dxFoiX) = 2tt' \ogK+^y- 4F3 1, 1, -, -; 2, 2, 2; 



3 3 



16 



^2 .^o. ^,^,2>2''-'-'-' ^2 



This agrees with (5.23) on the ABJM slice. To see this, one notices that 



a;W(z,^) 



2EE 

n=l k+l=n 



{2k + 21- 1)! 



y)2(UV 



{k\f{l 



2E 



4"(2n- l)!r(n+i 



n=l 



^T{n+lf 



21 ^n 



4z4i^3(l,l,^,^;2,2,2;16z 



is precisely the generalized hypergeometric function appearing in (5.31) 



(5.31) 



(5.32) 
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We are now ready to discuss the calculation of the planar free energy at strong coupling. We 
have, 

a^Fo(Ai,A2) = 2T:'^\ogK-TT'^J^\z,z). (5.33) 

After plugging the value of k in terms of A given by the series expansion (5.11), and integrating 
w.r.t. A, we obtain 

F„(A,B) = 14^A-V' + X:e-'^/, ('p] -'4{b-\)\ (5.34) 

where fi{x) is a polynomial in x of the form 

2£-3 

/,(x,/3)= V/f(/3)x^ t>2. (5.35) 



The coefficients /j. (/3) are Laurent polynomials in /3 of degree I, and symmetric under the 
exchange /3 -H- fi~^. We have, for the very first cases, 

/i(a;,/3) = -J(/3 + r'), 

/ (5.36) 

/2(x,/3) = ^(/32 + i6 + r') + f (/3 + r')'. 

The integration constant in going from (5.33) to (5.34) can be seen to be zero by comparing 
(5.34) with a numerical calculation of the integral (5.30) at intermediate coupling. 

The free energy in the planar approximation is given by rescaling (5.34) by the string coupling 
F = Qs'^Fq + 0{g^). This expression displays many interesting features. First, note that on the 
ABJM slice Ni = N2 the leading term 

_W2^2p/2 (537) 

displays the "anomalous" scaling N'^''^ in the number of degrees of freedom for a theory of M2 
branes, as was first derived from a supergravity calculation in [25]. The above calculation is a 
first principles derivation of this behaviour at strong coupling in the gauge theory. Usually, this 
behaviour is associated to the thermal free energy on M^, while (5.37) gives rather the free energy 
of the ABJM theory on S^ at strong coupling. However, a supergravity calculation of this free 
energy also leads to the A^^" behavior. We will show this now, and in particular we will match 
the numerical coefficient in (5.37).^ 

5.4 Calculation of the free energy in the gravity dual 



Consider type IIA theory on AdS4 x CP , and let us reduce it to the AdS4 factor as in for example 
[50]. The (Euclidean) AdS metric appropriate for a boundary theory on S"^ is 

ds^ = d/?^ + sinh^ pdf]2, (5.38) 

where dJl^ is the metric on an S'^ of unit radius. In this coordinate system, the boundary is 
at /9 — )■ 00. The free energy of the boundary CFT on S^ should be given, in the supergravity 



^We would like to thank Diego Hofman for very useful remarks on this calculation. 
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approximation, by minus the Euclidean gravitational action of the AdS„+i space — /AdS„+i- This 
action is given by a bulk term, a surface term, and a counterterm at the boundary [51, 52] 

-^AdS„+i = -^bulk + Isuri + -^ct, (5.39) 



with 



'bulk 



'surf 



1 

IQttGn 
1 

'SttGn 



d 



n+l. 



X 



dJ'x VhK, 



dX 



(5.40) 
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ct 



SttGn- 



(TxVh 



dX 



n 



1 + 



1 



2(n 



-11 + 



In these equations, Gn is Newton's constant, and R, K and IZ are the scalar curvature of the 
bulk, the extrinsic curvature of the boundary dX, and the scalar curvature of the induced metric 
h on dX, respectively. The counterterm action includes higher order corrections which are not 
relevant for the case of AdS4 and will not be considered here [52]. As our boundary dX, we will 
take the hypersurface p = pQ, and at the end of the calculation we must take po — )■ oo. The 
counterterms guarantee that the resulting action will be finite. 
The bulk action is easy to evaluate and gives 



-^bulk(PO 



SttGn 



vol(AdS4;/9o) 



where 



vol(AdS4; po) = vol(S^) / dp (sinh pY = 2tt' 

'0 



po 



1 3 2' 

— cosh(3/9o) - - cosh(/3o) + - 



(5.41) 



(5.42) 



It is easy to see that the surface term and the counterterms remove the divergences as po ^ oo, 
leaving only the term 47r^/3 in (5.42), and we find [52] 



Ihn /AdS4(Po) = 1^;^^- 
po->oo ^(-J^N 

If we now use the dictionary relating Newton's constant to the gauge theory data, 



2^/2 



Gn 



k'X^/', 



(5.43) 



(5.44) 



we find exactly the leading term in (5.34)! Of course, in order to obtain this result we have used 
the regularization provided by the counterterm integral in (5.40), and one could suspect that the 
matching depends very much on this regularization. However, this counterterm has been tested 
(or fixed) in an independent way in the calculations of [51, 52]. In particular, for n = 4 it leads 
to the matching of the Casimir energy of AA = 4 SYM on M x S^ , and for n = 3 it reproduces 
the standard mass of an AdS4-Schwarzschild black hole [51]. Therefore, the above calculation 
provides a genuine test of the AdS4/CFT3 correspondence. 

In Fig. 4 we show the exact result for the planar limit of dxFQ{X) in the case Ni = N2, as a 
function of A = N/k, and we compare it to the behavior of the supergravity prediction 



dxFoiX) w 27rV2(A - 1/24), A ^ 00. 



(5.45) 
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Figure 4: Comparison of tlie exact result for d\FQ{X) given in (5.29), plotted as a solid blue line, and 
the weakly coupled and strongly coupled results. In the figure on the left, the red dashed line is the 
supergravity result (5.45), while in the figure on the right, the black dashed line is the Gaussian result 
(5.46). 

We see that the strong coupling behavior gets triggered for values of the coupling A ~ 0.2. For 
A — )• 0, the behavior of the prepotential is dominated by the Gaussian, weakly coupled result 

(4.4) 

dxFoiX) « -Svr^A flog ('^] - l] , A ^ 0. (5.46) 



A second aspect to notice is that the supergravity result (5.34) has corrections which are 
exponentially suppressed. The exponential is of the form 

.1) 



where 



is the area of the CP^ two-cycle in 

tions multiplies (at each order in A""*^") the polynomial /^ (/3) in P,/3^^ 

contributions schematically of the form 



(5.47) 



27rV2A (5.48) 

Also, notice that each of these exponential correc- 

Therefore, we have 






-n+(A(CPl)+27riB)-n_(yl(CPl)-27riB) 



(5.49) 



This is precisely what one should expect for a gas of n^ instantons and n_ anti-instantons 
in a 0" model on CP^, where the (anti)instantons wrap the CP^ cycle. Notice that this kind of 
corrections are made possible by the non-trivial topology of two cycles in CP^, i.e., by the fact that 
62(CP^) = 1, and as such they are absent in AdSs x S^. Some aspects of these string instantons 
have been studied in [53] . It would be interesting to test in detail the possible connection between 
these string instantons and the exponentially suppressed corrections to the planar free energy. 

These instanton corrections are also present in the Wilson loop result (5.16), again with an 
infinite series of corrections. This can be compared with the case of AA = 4 SYM in 4d, where the 
asymptotic large coupling expansion of the Gaussian matrix model (2.32) has a single instanton 
correction which can be explicitly identified with a second saddle point solution in AdSs x S^ 
[48, 54]. 

Finally, we note that when Ni ^ N2, the planar free energy (5.34) includes an imaginary 
term proportional to (B — 1/2)'^, which is derived by the weak coupling calculation (4.5). In CS 
theory such a term is related to framing [12]. It would be very interesting to derive this phase 
in type IIA string theory. 
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6. Conifold expansion 

The expansion around the conifold locus corresponds to a region in the moduli space of the 
ABJM model where one of the gauge groups has finite coupling, while the other one is weakly 
coupled. In the lens space matrix model this corresponds to one 't Hooft parameter being small, 
and the other of order 1. In this section we will study this regime from three different points 
of view: the exact planar solution in terms of periods and Picard-Fuchs equations, the matrix 
model, and the gauge theory. 

6.1 Expansion from the exact planar solution 

We can use the exact planar solution to calculate various physical quantities near the conifold 
locus. For concreteness, we will expand around t2 = but with ti arbitrary. The first ingredient 
we need is an expansion of the global coordinates of moduli space. It turns out that the most 
convenient method is based on the expressions for the periods (2.14). The locus where ^2 = is 
the line 

C = 2/3-2, (6.1) 

where the cut (—6, —1/b) collapses to the point Z = —1. The derivative of t2 w.r.t. ( can then 
be computed in terms of residues at this point by expanding the expression in (2.16): 

-^ = 1^^J^^ (Z + l)2m ' (6-2) 

where Hk{Z^ f3) are regular at Z = —1. This gives a series for ^2 in powers of C — 2/3 + 2, 

which can be easily inverted to 

C = 2p-2-Ay^t2 + l{l-^)4 + ^-±^^^4 + Oi4). (6.4) 



As a nice application of this expansion, we can compute the VEV of the 1/2 BPS Wilson loop 
around the conifold point, which is given in (2.35). Using the dictionary (2.23), (2.24), we find 

e-'''^{W'^^^)g=o = 2 sin(7rAi) + 27rA2 (2 - cos(^Ai)) + tt^X^ sin(7rAi) 

1 (6.5) 

+ -vr^A^ (l - 5cos(7rAi) + 3cos^(7rAi)) + 0{X^). 

As A2 — ;• 0, we recover the result for a Wilson loop VEV in U{Ni) CS theory. In the conifold 
expansion we are then regarding the ABJM theory as a perturbation of U{Ni) CS theory at 
strong coupling. 

The above result can be also obtained by solving the Picard-Fuchs equation around a point 
in the conifold locus. Let us choose for example the symmetric conifold point (3.32), with B = 1 
and K = 4. This corresponds to the point in the conifold locus with 

Ai = ^, A2 = 0. (6.6) 
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The appropriate global coordinates near this point are (3.33). We find that A2 is a period solving 
the PF system (3.34) and with leading behavior 

X2 = -^{y2 + yi/2) + 0{y^). (6.7) 

One finds the expansion 

X2 = -{B -if-'^iB- If + f ^ - ^{B - If + ^{B - If) (K - 4) 
4^ ^ 96^ ' \8tt 32^ ^ 1536^ V 

+ ( — + -^{B - if - —{B - if) {k - 4)2 + OUB - if) + OUk - 4)3), 

V 1287r 1024^ ^ 8192^ V ^ ^ VV JJ^W JJ^ 

(6.8) 

which is inverted to 



-2 



1\' vr^ /. 1 



k=A-2-kU Xi- -] +— Ai-- +7rA 




2J 6 V" 2J ' ^y ' "V" 2J 3 V" 2J I (6.9) 

+ 0(Ai)+0((Ai- 1/2)5). 

This is indeed the expansion around Ai = 1/2 of (twice) the series in the r.h.s. of (6.5). 

Once we know the expansion of the global coordinates, we can consider other quantities in 
the model, like the genus g free energies. The conifold expansion of Fg{ti, ^2) has the form 

F,(Ai,A2) = F«(A2) + J;fW(Ai)A^, (6.10) 

n>0 

where F5'(A2) is the free energy of the U{N2) Gaussian matrix model, and each coefficient 
Fg {Xf can be obtained as an exact function of Ai. Of course, 

Ff(Ai) = Ff(Ai) (6.11) 

is the genus g free energy of the CS theory on S^. When </ = 0, the expansion (6.10) can be 
computed from the exact planar solution in various ways. One can for example use the Yukawa 
couplings (3.37) expanded around the conifold locus in order to compute the third derivatives 
of Fq, or use the modularity properties of the solution discussed in [41, 31]. In any case, for the 
first few functions one finds the following results: 

Fo(')(Ai) = 27ri U^Xl + 2Li2( - e'^^i) - 2Li2( - e""^^!^ 



F^^\xf = -27r3iAi +8^2log (^cos (^)) , (6.12) 

Fri^i) = ^- + y (3cos(7rAi) - 5) tan (^^j . 

6.2 Conifold expansion from the matrix model 

It is easy to implement the conifold expansion directly in the lens space matrix model. To do 
that, we notice that it can be written as two interacting Chern-Simons matrix models on S^. We 
recall that the CS matrix model on S^ , first considered in [9] , is defined by the partition function 

z..(iv,*)=^/nt^n(2--''('^))'''"-=-'''- «5-") 
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This is a one-cut matrix model [21]. It can be obtained from the lens space matrix model when 
one of the two cuts collapses to zero size. In the Z plane the endpoints of the cut are given by o 
and a~^, where 

a = 2e*- l-2e*/Ve*-l. (6.14) 

Let us consider the following operator in this model: 

y^{vi) = 2 J] log (2 cosh (^^^)) • (6.15) 

The lens space partition function (2.3) can be calculated in two steps. In the first step, we 
compute 

Zi(z.,) = (e^^'^^))^^ (6.16) 

where the subindex A^i indicates that this is an unnormalized VEV in the S^ CS matrix model 
with gauge group U{Ni). In a second step, we calculate 

^L(2,l) = {Zlil^j))N2 (6.17) 

in the CS matrix model with gauge group U{N2)- To obtain the conifold expansion, we calculate 
Zi{vj) and we expand it in Qs and around Vj = 0. Each term in this expansion can be computed 
exactly as a function of the Kahler parameter ti, since the CS matrix model can be solved exactly 
in the 1/iV expansion. The resulting double series in Qs and Vj is then regarded as an operator 
in the CS matrix model with group U{N2), which we expand around the Gaussian point as in 
[9, 10], i.e., we expand the sinh measure around Uj = 0. The partition function .^^,(2,1) is then 
computed as a VEV in the Gaussian matrix model. This procedure gives a method to compute 
the expansion (6.10) directly in the matrix model. 

To illustrate this procedure, let us calculate -^0(^1)^2) at first order in t2. In this computation 
we will denote 

Ui = diag(e^')> U2 = diag(e''0. (6.18) 

The expansion around z^j = of the operator W^i^j) reads 

Ni N2 Ni 

W(i.,) = 2iV2X;iog[2cosh(|)]-j;z.,J]tanh(|)+0(^|). (6.19) 

1=1 j=l i=l 

The average of the second term in the U{N2) matrix model vanishes (since it is odd in I'j), while 
higher order terms are at least of order t^. The first term can be written as 

2 ^ log [2 cosh (^)] =2TTlog{l + Ui)-^fii. (6.20) 

j=l i=l 

Therefore, in the planar limit and neglecting terms which contribute at order t|, we have 

logZ,{u,) « -^ (Tr log(l + t/i))jv, (6.21) 

9s 

since the second term in (6.20) is odd in fii and its VEV vanishes. We then find, 

Fo{h,t2) = Ff (ti) + 2t2gs (Tr log(l + Ui)) + 0(4). (6.22) 
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The VEV in (6.22), which is now norniahzed, can be computed in terms of the resolvent of the 
CS matrix model, and similar computations appear in [55, 14] in the context of large A^ instanton 
corrections. In fact, it follows from (8.28) and (8.30) that the VEV in (6.22) is given by —g{—l), 
where g{Y) is computed in (B.2). The final result for the linear correction in ^2 is 

^ + I + Li2(e-*^) - 2Li2(e-*i/2) + 2U2{-e''^/^). (6.23) 

Using dilogarithm identities, this agrees with -^ig (^^i) ™ (6.12). It is interesting to point out 
that, in the context of CS theory on the lens space L(2, 1), this function is essentially the action 
of the large N instanton corresponding to the flat connection 

U{N) -^ U{Ni) X U{N2), N2 < A^i, (6.24) 

as shown in [14]. In the matrix model, this action is obtained by tunneling N2 eigenvalues from 
the first cut to the second cut. 

We can also calculate the conifold expansion for the VEV of 1/6 and 1/2 BPS Wilson loops 
directly in the matrix model. We want to compute 

{w'J')=gs{TYUi)n2,i)- (6.25) 

We will again perform this computation in the planar approximation and at linear order in t2- 
At this order we can compute instead the normalized average of the operator 

^ /VV(.)\ ^^^ =(Tr^i) + (Tr^iW(i^,))^^^ + --- (6-26) 

in a Gaussian matrix model for the Uj. In the last line, all VEVs are normalized VEVs in the S^ 
CS matrix model. By completing the square of the Gaussian weight we derive 
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TVt/i(^j;/.j\ = |(Trt/ie^-S">' 



= gs{TvUi). (6.27) 

j=0 



We then find, at this order, 

(I^;/')5=o = 9s{TvUi) + t2 (2{TiUi Trlog(l + Ui))^'^ - ffs(Tr iJi)) + 0(4). (6.28) 

The connected correlator 

°° (-^V 
(Tr C/i TVlog(l + C/i))(=) = -^ ^-^(Tr f/i Tr t/f)(") (6.29) 

e=i 

can be computed by considering the (partially) integrated two-point function (see for example 
[56]) 

fdpWo{p,q) = -J2—^,{TrU^^TTUr)^^^ (6.30) 

and extracting the coefficient of q'~'^. We have, 

dpWoiP,<l) = , ' , (l - J ^fl^^~"']l ] + , ' (6.31) 

"^^'''^ 2{p-q)\ ^ {q-a){q-a-^) 2^{q - a){q - a^^)' ^ ' 
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which includes the appropriate integration constant. We find, after changing p — t- —p, 

-E^^Tr;7iTrC/f)W = - (^a + a'' + 2p - 2y'ip + a)ip + a-^)) . (6.32) 

e=i ^ 

When p = 1 this gives 

- J2 ^T^(Tr Ui Tr C/f )(^) = e*i - e^^/^^ (6.33) 

e=i ^ 

Notice that this is an infinite sum of correlators in the CS matrix model. Since 

{TtUi) = ^^^, (6.34) 

gs 



we finally obtain, 



{w'J')g=o = e'^-l + t2(e'^/'-l) +0{4) 



^*i/2 I 2 sinh ^ + t2 ( -2 + 2 cosh ^ ) + 0{tf 



(6.35) 



Since this is a Wilson loop only in the first group, the framing prefactor depends only on the 
first 't Hooft coupling. 

The 1/2 BPS Wilson loop is obtained by subtracting 

(Tr f/2)L(2,i) = N2 + 0{4) = ^ + 0{tl). (6.36) 

9s 

We find, 

^-(..-ft.)/2^^V2^^^^ = 2 sinh (I) + t2 (-2 + cosh (|) ) + 0{tl). (6.37) 

This is the result (6.5) obtained from the conifold expansion after using the dictionary (2.22). 

6.3 On the near Chern— Simons expansion of ABJM theory 

In the matrix model the conifold locus corresponds to vanishing of one of the two cuts, where 
the lens space matrix model can be written as a perturbation around the matrix model for 
Chern-Simons on S^. Here we want to explore this limit in the original 3-dimensional theory. 

In the strict limit we have the theory with A'^2 = and Ni ^ 1 and arbitrary Ni/k. 
In this limit all the fields charged under the second gauge group, i.e., its gluons and all the 
bi-fundamental fermions and scalars are removed. Consequently, ABJM theory simplifies dra- 
matically and reduces to topological U{N) CS. The only observables in the theory in this strict 
limit are Wilson loops, and they are given by the standard CS answer, which is exact in Ai (and 
l/Ai). 

One can try to perform a systematic expansion around this point in a perturbative expansion 
in A2. One keeps A2 <C Ai, but if desired, can still assume the planar approximation, ignoring 
also the I/N2 corrections. 

It is convenient to draw the Feynman graphs in double-line double-color notation, one color 
for each group. At the first non-trivial order in A2, only graphs with a single index loop of U{N2) 
are included. An arbitrary number of gluons of U{Ni) are allowed. Let us propose the following 
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calculation procedure: First ignore all U{Ni) gluons and enumerate all remaining graphs. They 
are a very restricted subset, which can be identified very easily. 

In order to dress them up with the U{Ni) gluons we write the proparators for the bi- 
fundamental fields as a path integral over all trajectories in space. As charged object, these 
paths will effectively be Wilson loop in U{Ni), which can be calculated exactly in CS theory. 
Since this theory is topological, the result of adding all the gluons, is independent of the path of 
the bi-fundamental fields. One can then do the usual path integral for these fields and find the 
regular scalar and fermion propagators. 

The statement in the previous paragraph fails in a subtle way. The correlation function of 
Wilson loop operators in CS theory does not depend on their geometry only as long as their 
topology — the knotting and linking numbers — are kept fixed. Therefore one has to modify the 
above statement, and sum over all possible topologies of the paths of the bi-fundamental fields 
accompanied by the relevant knotted/linked Wilson loops. Unfortunately, we do not have an 
a-priori method of determining the weight that should be assigned to the different topologies.^ 

As an illustration, let us consider the 1/6 BPS Wilson loop (whose Feynman rules are simpler 
than the 1/2 BPS one) and examine its perturbative expansion about the conifold locus. The 
Wilson loop is given in our normalization by [15, 16, 17] 

W^^^ = Qs TrPexp / ( iA^x^ + '^\x\M^jCicA ds. (6.38) 

x^ parameterizes a circle in M^ (or S'^), A^ are the U{Ni) gluons, C/ and C^ are the bi- 
fundamental scalars and Mj = diag(l, 1, —1, —1) is a matrix in flavor space, which is required to 
make this object BPS. 

At order 0{X^), this is simply a Wilson loop of CS, whose planar expectation value (ignoring 
framing) is 

{Wll^)g=0 = 2isin7rAi + ©(As). (6.39) 
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Figure 5: Several Feynman graphs which may contribute at order A2 to the 1/6 BPS Wilson loop with 
gluons stripped. The big circle is the Wilson loop, dashed lines are bosons and the solid line a fermion, 
all presented in double-line, double-color notation. 

After stripping away the gluon lines, there are still an infinite number of graphs involving 
bi-fundamental fields. Examples are shown in Fig. 5. In the example drawn, there is a single 
scalar or fermion loop. The scalar loop can "touch" the Wilson loop at an arbitrary number of 
points, due to the scalar bilinear term in (6.38). There are extra graphs which are not drawn, 
with fermionic tadpoles on the scalar lines, or vice-versa. 



* Wilson loops arise out of dressing propagators of matter fields also in [57] . In that case the path is fixed to a 
collection of light-like segments, due to the singularity in the Minkowski-space propagator. 



32 



By explicit calculations [15, 16, 17], all the connected graphs illustrated (Fig. 5c, 5d, 5e) 
vanish in dimensional regularization. The same can be argued for higher order graphs of this 
form. Likewise, one would not expect tadpoles to contribute. We are left therefore with the 
first two disconnected graphs, which become connected once gluon lines are added. Indeed, the 
only non-vanishing graph that was thus far calculated is the one-loop correction to the gluon 
propagator (Fig. 5a, 56 with two extra gluons), which accounts for the 0(A2Af) term (which 
with our normalization is 2-loops) in the explicit answer (4.3).^ 

We can compare this to the explicit calculation in the matrix model above. The essential 
part of the expression for the Wilson loop at order 0{X2) (6.28) is the connected correlator of 
two Wilson loops. One of them is the original Wilson loop and the other came from expanding 
the cosh term in the matrix model (6.20), which arises from integrating out the bi-fundamental 
matter. So this agrees with the identification of the contribution as coming from the bubble 
graphs. Moreover, what we see in the matrix model is that one should sum over multi-winding of 
this second Wilson loop, with a weight l/l. This corresponds in the physical theory to summing 
over all possible topologies for the scalar and fermion bubble. As mentioned, we do not know 
how to derive this factor of l/l from perturbation theory, but it is given to us by the explicit 
matrix model calculation. 

It was noted in [26] that in this limit of ABJM theory the spectrum of local operators also 
simplifies and the spin-chain hamiltonian becomes short-range. A compelling conjecture for the 
mysterious function h{X) in that limit was also presented there. It would be interesting to explore 
this limit further and learn how to do this sum over topologies for other observables. 

7. Modular properties and the genus expansion 

In this section we provide an efficient, recursive method to compute the 1/A^ corrections to the 
free energy in the case Ni = N2 = N. This is based on the modular properties of the solution 
and the technique of direct integration of the holomorphic anomaly equations. The method 
determines a priori the full 1/A^ expansion. In practice it is quite efficient and it makes possible 
to calculate the Fg corrections for high genera. This is then used to estimate non-perturbative 
effects in the large A^ expansion. 

As noted in [41], we can use the relation between the local Fq theory and Seiberg-Witten 
theory to write all the quantities in the model in terms of modular forms. This representation 
becomes particularly useful when we restrict ourselves to a one-parameter model, as it was shown 
in a different context in [58]. When Ni = N2, 13 = 1 and the modulus u becomes simply 

« = 1 + -. (7.1) 

In Seiberg-Witten theory, u is related to the modular parameter r of the Seiberg-Witten curve 
by 



a- '^ ^ "^ (t) = 1 - 32g^/^ + 256g + • • • (7.2) 

where q = e^'^'^. This formula can be inverted to 

r-i^^ (73) 



®This graph has a divergence that can be removed by including the double scalar exchange graph (Fig. 5c) . In 
dimensional regularization the finite part comes only from the gluon graph. 
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therefore we see that the modular parameter r is related to the specific heat of the theory through 
(5.26). Let us now introduce the quantity 

This is proportional to the third derivative of the genus zero free energy, therefore to the Yukawa 
coupling Caaa- More precisely, we have 

dlF^iX) = -S^ht (7.5) 

Therefore, the planar content of the theory can be elegantly encoded in terms of modular forms 
on the Seiberg-Witten curve. 

One powerful application of the modular properties of the AB JM theory is the determination 
of the higher genus corrections to the free energy, Fg{\). These can be obtained in principle from 
the matrix model (2.1), or equivalently from the formalism of [59] (appropriately modified as in 
[60, 61]). However, as emphasized in for example [28, 62, 58], this formalism is not very convenient 
to do calculations at higher genus. One should rather use the fact that the Fg are quasi-modular 
forms that can be promoted to non-holomorphic modular forms. The resulting non-holomorphic 
objects satisfy the holomorphic anomaly equations of [30], as shown in [28, 63], and these can be 
in turn solved with the technique of direct integration developed in [28, 29, 31, 58] for local CY 
manifolds and matrix models. 

The basic strategy of direct integration is the following. First, we assume an ansatz for Fg 
of the form 

F,{t) = e'~'fa{r) (7.6) 

where 

39-3 

fair) = J2 Elir)ci'\r) , g>2, (7.7) 

fe=o 

is an almost modular form of weight 6g — 6, with respect to a monodromy group T C SL{2,Z). 
Fglr) can be promoted to a non-holomorphic modular form Fg{T,f) by changing 

E2{t) ^ E2iT,f) = E2{t) - — ^. (7.8) 

The resulting Fg (r, r) satisfies the holomorphic anomaly equations of [30] , which govern their 
anti-holomorphic dependence. Since this dependence is contained in E2{T,f), these equations 
govern the E2 content of Fg. This means that the coefficients c^ (t), which are modular forms 
of weight 6(7 — 6 — 2k, can be obtained recursively for A; > if one knows the lower Fg. In order 
to write down the recursive equation, it is useful to introduce a covariant derivative dg taking a 
form of weight k into a form of weight k + 2: 

^^ = 9r + l^ (7.9) 

Then, the holomorphic anomaly equations lead to 

^ = -3 Uf9-i + 3 ^ d^/9-1 + Ed^/'-d^/^-} ' 9>2. (7.10) 
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If Fgi are known, with g' < g, the above equation determines all the coefficients cf (t) in fg, 

with the exception of Cq (r), which plays the role of an integration constant. This coefficient is 
a holomorphic form of weight 6g — 6 and it is called the holomorphic ambiguity. 

In order to fix the holomorphic ambiguity we need two pieces of information. The first one 
concerns its functional dependence. Since Cq (r) is a modular form w.r.t. some monodromy 
subgroup, it belongs to a finitely generated ring. This means that it is determined by a finite 
number of coefficients, which typically grows with g. The second piece of information comes from 
boundary conditions at singular points in moduli space. A very powerful boundary condition for 
matrix models and local Calabi-Yau manifolds is the so-called gap condition, discovered in [28] 
and further used in [31, 58] to fix the holomorphic ambiguity. According to the gap condition, 
near certain points pi in moduli space, parametrized by a flat coordinate tj, the genus g free 
energy behaves as 

^i^^ = A + 0(1)- (7-11) 

The superscript (i) means that the genus g free energy has to be transformed to the duality 
frame which is appropriate for the i-th singularity, as it is well-known in special geometry. The 
"gap" refers to the absence of singular terms t~ with < k < 2g — 2 in the local expansion near 
ti = 0. The vanishing of these terms provides boundary conditions for Cq (r), and in some cases 
it fixes them completely. 

In our case, the relevant ring is that of r2 modular forms which is generated by the theta 
functions 

b = 4{T), c = ^I{t), d = ^i{T). (7.12) 

Since c = b + d, only two of them are independent, and we will choose b and d. Using standard 
formulae in the theory of modular forms, one finds 

(7.13) 











dri b - E2 










i 4 


as well as 












d^b = 


b^ + bd 
3 ' 


d^ibd) -- 


{bd)b 

= 6 ' ^«^ 



-^{-El + 2bE2-E^). (7.14) 

The modular expression for the genus one free energy is known [41] and reads 



therefore we have 



Fi = - 


-log?7(r). 


d^/i^ 


E2 

24' 



(7.15) 



(7.16) 



These are all the ingredients needed for the recursion. The holomorphic ambiguity can be written 
as 

43)(r) = ^ afVd^s-^-^ (7.17) 

and it involves 3(7 — 2 unknowns. Let us see how we can fix these by looking at the behavior near 
the three singular points of moduli space. 
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At the orbifold point, the Fg are the genus g amphtudes of the super- matrix model (2.1) 
with A^i = N2- Their leading behavior near A = is governed by two copies of the Gaussian 
matrix model, therefore they behave as 

^i°^(A) = — ^(27riA)2-2^ + 0(1). (7.18) 

" g{^9 - 2) 

This gives g — 1 conditions, since the ansatz (7.17) for the holomorphic ambiguity only involves 
even powers of A. 

The symmetric conifold point zi = Z2 = 1/16 is related to the orbifold point through an 
S'-duality transformation. The appropriate global coordinates near this point are given in (3.33). 
In the ABJM slice one has 

yi = 0, y2 = y = l-^. (7.19) 

lb 

The following period is a good local, flat coordinate near the symmetric conifold point: 

oo 

t = y , \ , ;/"+S (7.20) 

where 

"" ~ 'an\ ^^\k J \2kJ V n-k }\2n-2k, 

\n) fc=0 



fn + l)24™' 

n=0 



Y^ f2k\ fAk\ (2n - 2k\ /4n - 4A; 



(7.21) 



F(=)(t) = ^ ,fr ^, (^ ' + 0(1). (7.22) 



It was noticed in [31] that the genus g amplitude at the conifold point behaves like 

B2g (t 

2g{2g-2) \2\ 

This fixes 2g — 2 conditions. 

Finally, the large radius point is related to the orbifold point by an STS transformation. 
The genus g free energy is the generating function of Gromov-Witten invariants of the local Fq 
geometry in the slice Ti = T2 = T. More precisely, one has 

d>l 

where 

Nd,s= J2 ^'^i-'^^'S (7-24) 

di+d2=d 

is a sum of Gromov-Witten invariants at genus g, iV^^^dj.g; of local Fq (the degrees di, d2 corre- 
spond to the two Kahler classes of this geometry). Since (7.23) is a power series in Q with no 
constant term, we obtain one extra condition, which, together with the g — 1 conditions from 
the orbifold point and the 2g — 2 conditions of the conifold point, completely fixes the 3g — 2 
unknowns in the holomorphic ambiguity. 

Let us see how this works in some detail when g = 2. The integration of the holomorphic 
anomaly equation gives. 



/2 = ^ T^ ( -^ ^2 + ^bEi - 2E,E2) + 4'V), (7.25) 



1 W 5 

3 " 242 V 3^^ ' ""^^ -^'i^y 
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(2) 

where Cq (t) is of the form (7.17). The expansion around the orbifold, conifold, and large radius 
points read, respectively, 

^(c) ^ 5 + 1296af -l-864(12af + 15af) ^ (726) 

2 ^ '^ 1296*2 10368t ^ ^' 

Fr\Q) = -4^ (3 + 1728(«f - af) + af - af )) + 0(Q). 
Imposing the conditions (7.18), (7.22) and (7.23) we fix 



«(")-— ^ a(2)-^^ a^')--^ a^"^-— f7 27l 

"0 " 25920' 1 " 17280' ^ " 3455' «3 " 3240" y'-^'> 



We finally obtain 



,(0) _ 1 f 5^3^,,rf .^^A . 16b^ + 15db^ + 2l£b + 2d^ 



Fk = r^ — E^ + 3bE^ - 2E4E2 + ^ ^^ ^ . (7.28) 

Since r depends on A through (7.3) and (5.5), this gives the exact expression for the genus two 
free energy on S^ in the ABJM model, for any value of the 't Hooft coupling. 

Notice that the modular ring appearing here and parametrizing the holomorphic ambiguity 
is different from the one appearing in Seiberg-Witten theory [28, 29] or in the cubic matrix 
model [58]. This is due to the fact that, although the curves are the same, the meromorphic 
forms defining the theory are different. 

Using this method, we have computed the free energies up to high genus. The strong coupling 
behavior of Eg is of the form 

FJ°)(A) ~ -X'2-9^ A ^ 00, g>0. (7.29) 

We have also used these results in order to investigate the large order behavior of the 1/A^ 
expansion. We have found that 

F,(A) = i-iy-' (fJ°) - -^-^(2vriA)2-^^) (7.30) 

behaves at large g as 

F,(A) ~ (25)! |^(A)|-'^ cos {2g0{X) + 5(A)) . (7.31) 

In this equation, the angle ^(A) satisfies ^(0) = tt/2 and ^(A) 7^ for all A > 0, while (5(A) is a 
function of A (see for example section 2 of [64] for more details on the large order behavior of 
the genus expansion). The sign (—1)^"^ is included in (7.30) since in the physical ABJM theory 
the coupling gg is imaginary. The large order behavior (7.31) indicates that the singularities of 
the Borel transform of Fg{X) which are closest to the origin are located at ib^(A), where 

A(A) = |.4(A)|e'^(^). (7.32) 

Since ^(A) does not vanish, none of them lies on the positive real axis. This strongly suggests 
that the 1/A^ expansion of the free energy is Borel summable for any A > 0. 
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The large order behavior of the genus expansion (7.31) is similar to the one found for Chern- 
Simons theory on S^ in [65], and it should be governed by a large N instanton with action ^(A). 
It would be very interesting to identify this instanton and compute A{X) analytically, both in 
the gauge theory and in the string theory dual. The factorial growth, found here by explicit 
calculation in the matrix model, agrees with the expected behavior for the genus expansion in 
string theory [32]. 



8. More exact results on Wilson loops 

In this section we elaborate on the results of [7, 8] and we obtain more exact results on Wilson 
loops. 

8.1 1/A^ corrections 

The higher genus corrections to the VEV of 1/2 and 1/6 BPS Wilson loops can be computed in 
terms of the higher genus corrections to the resolvent of the matrix model. The resolvent has a 
genus expansion of the form 

oo 

u;{z) = Y,9?^9i^)- (8.1) 

3=0 

In the same way, the density of eigenvalues has a large N expansion of the form 

oo 

p{z) = Y, 9?Pg{z) , p{z) = p(i) (.) + p(^) {z) . (8.2) 

9=0 

The Pg (z) (with i = 1, 2) have their support on the intervals Ci, and they can be obtained by 
the discontinuity of ujg at the cuts as in (2.28). 

The genus expansion of the expectation value of the 1/6 BPS and 1/2 BPS Wilson loops 
follows the expressions in (2.30) and (2.31) with the appropriate term in the expansion of p^^'(Z) 
and uj{Z). 

The first step is therefore to compute ojg{p). This calculation can be done with the recursive 
techniques developed in the matrix model literature starting with [56] and culminating with [59]. 
We will perform an explicit computation for g = 1- Calculations for g >2 are in principle doable, 
but they become complicated. 

A convenient formula for uji{p) for an algebraic resolvent was found in [66]. To write this 
formula, we write the discontinuity of the resolvent (also called spectral curve in the matrix model 
literature) as 

y{p) = M{p)y/a{p), a{p) = (p- xi){p- X2){p-X3){p-X4). (8.3) 

M{p) is sometimes called the moment function. Then, one has 

^^iv)-^^t{j^_.^^.c), ,8.) 
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where 



Ai 
Bi 

Ci 



1 1 

16 M(xi)' 

1 M'{x,) 
16 M2(xi) 

1 1 



+ 8M(x,)r"^ E^^_^J' 



^E 



a,- — Oi 1 M'(xi) Oi 



(8.5) 



48 M(xi) V ^7 - Xi 16 M2(x 



and the ai are given by 



1 



ai 



02 



as 



a4 



(xi - X2) 
1 

(X2 - Xi) 

1 
(X3 - X4) 

1 



8M(x,) 

_ (X4 - X2) -E(fc) 

(x4 — xi) K{k) 

(X3 - Xi) E{k) 
(X3 - X2) Kik) 

_ (x4 - X2) E{k) 
(xs - X2) K{k) 

^ (X3 - xi) E{k) 
(x4 — xi) K{k) 






(8.6) 



(X4 - X3) 

where the modulus of the eUiptic functions is 

^2 _ (a^i -a:2)(x3 - X4) 
(Xi -X3)(X2 -X4)' 



k' 



(8.7) 



These expressions differ from the ones in [66] in a permutation of the roots, as explained in [67]. 
The overall factor of 4 in (8.4) is due to the fact that our resolvent has a different normalization 
than the one in [66]. 

Although the resolvent of the lens space matrix model (2.10) is not algebraic, its discontinuity 
can be written in the form (8.3) with 



and 



ct(p) = f{pf - 4/?V, f{p) =p^-Cp+i 
2 



M{p) 



■ tanh' 



-1 \A(p) 



(8.9) 



This form of the spectral curve is typical of the mirrors of toric geometries [60, 61]. The branch 

points are 

11 , , 

xi = -6, X2 = -7, X3 = -, X4 = a. (8.10) 

b a 

Using these expressions, it is possible to compute the integral 



«'°>'- 4„ 



— / uji{Z)ZdZ 



i.n) 



in closed form, in terms of elliptic functions E, K and the elliptic integral of the third kind 
n(n, /c), with 



n 



(g^ - 1)6 
(1 + ah)a 



^.12) 



39 



One finds the rather comphcated expression 



(</').=! 



1 



127rV^63/2(i + ab){a'^ - 1)2(62 _ i)K 



-3{b-2a + a^b){l + abyE^ + 



a(l + a 



b + 0^(4 + Aa^ - a^)b - 4a(l - 3a=^ + a^)b^ - a^{l + 0^)6^(2 + b^) + a(l - Sa^ + 0^)6^ 



+ (6^(1 + 6a^ + a^) + 4a(l + a2)(62 - 1) + 6(3 - Ua^ + 3a^)) (1 + abfE K 
(a6-l)(a2-62) 



K' 



+ 



12^(a6)3/2(l + a6)A:4i^2L 



6^;^ + 4(2 - k^)EK - (2 - 2^ + r) iiT' 



n. 



(8.13) 

To check this formula, we expand it around the weakly coupled point Ai = A2 = 0. After using 
the inverse mirror map given by (4.1) we find 



{W'j\=i 



vri 



vr 



TT 



vr 1 



^A, + -Xi + -A,A, + -A^ + -Af A, 



vr^i 



ttH 



-Aia: 



1A2 



|.Af + |^A?A, + ^A?A^-^A,Ai + 0(A^) 



^.14) 



We can test this expansion with a perturbative calculation in the ABJM matrix model. At order 
0{gg) we have found, 

(</') = !- 



27riAi 
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5760 ^ 192 32 ^ 5760 



(8.15) 
It is straightforward to see that this agrees with (8.14). 

The 1/A^ correction to the 1/2 BPS Wilson loop is much easier to obtain, since it can be 
computed as a residue at infinity. We have that 

4 



4 

Z2 



U^l{Z) = ^Y.^^ + 0{Z-^), 



^.16) 



where the Cj are given in (8.5). We find, at weak coupling 

^2 



{W'j'),=i 



TTl 

12 



vr^i 



(Ai + A2) + ^(A? - A^) + ^(A? + A^^ 



12 



18 



57r'^i 



AiA2(Ai — A2) 



|^(At-At) + ^AiA2(A?-Ai) + 0(A5). 



^.17) 



At strong coupling we find (we consider for simplicity the ABJM slice) 

1 3 + 2 log^ K - 4 log K 



{<%=! 



-K + 0(1) 



(8.18) 



24i log2 K 

The leading exponent is exactly as at genus zero (5.16), representing the same minimal surface 
with an extra degenerate handle attached. Its effect is to modify the one- loop determinant, which 
(with our normalization and ignoring instantons) can be written as 



{W, 



l/2v 

D '9=1 
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12 67r\/2A 167r2A 



,7rV2A 



A — )■ 00. 



(8.19) 
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8.2 Giant Wilson loops 

It has been argued in [33, 68, 34, 35] that a D-brane probe in AdSs x S^ represents an insertion 
of a Wilson loop in the dual 4d A^ = 4 SYM with a large symmetric or antisymmetric represen- 
tation (in the case of D3 branes and D5 branes, respectively). These "giant Wilson loops" are 
characterized by a representation with n boxes, and one considers the limit 

Tl 

n, N ^oo, — fixed. (8.20) 

In terms of the Gaussian matrix model of the Wilson loops in that theory, the giant Wilson loop 
in the symmetric representation is represented by an additional eigenvalue outside the cut and 
the antisymmetric representation by a "hole" in the original cut. 

Let us review now the known D-brane solutions which could be relevant for ABJM theory. 
The usual 1/2 BPS Wilson loop in the fundamental representation is described by a string with 
world- volume AdS2 C AdS4. In M-theory it is an M2-brane wrapping also the orbifold cycle 
on S^/Zfc. When considering k/2 coincident M2-branes (or k, when it is odd) the M2-brane 
solution develops an extra branch, where the circle becomes a linear combination of the orbifold 
direction and a contractible circle in AdS4 [69]. In type IIA these configurations are D2-branes 
with world-volume AdS2 x S^ C AdS4, where the radius of the S^ is a free modulus. From the 
M-theory point of view these are continuous deformations of the system of k/2 coincident M2- 
branes describing a Wilson loop in a k/2 dimensional representation. In the field theory they are 
the vortex loop operators of [18], which have a description as semi-classical field configurations 
and carry the same charge as k/2 Wilson loops. 

These solutions have further moduli associated to rotations away from the orbifold cycle 
inside S'^/Z^. Such M2-brane configurations preserve 8 supercharges (1/3 BPS) [15, 18]. 

There is also a known family of D6-brane solutions which were argued in [15] to represent 
the 1/6 BPS Wilson loops in anti-symmetric representations. The action for this D-branes is (for 
iVi = N2) 

5B6 = -vr^/^"^^, (8.21) 

which matches that of n strings for small n and has the n ^ N — n symmetry of the antisymmetric 
representation. In the matrix model these D6-branes should correspond to creating a "hole" in 
one of the two cuts, splitting it in two. 

We turn now to the lens space matrix model and try to find the appropriate description 
for these objects, and in particular the 1/2 BPS vortex loop operators. As pointed out in [70], 
the calculation of Wilson loops in the matrix model in this limit can be done in a saddle-point 
approximation. We will now reformulate the arguments of [70] and adapt them to the lens space 
matrix model. 

We will focus on the case of 1/2 BPS Wilson loops, where we want to calculate 

W:i = {TrT^vU), r/ = ±l, (8.22) 

where U is the same matrix as in (2.9) and 7^„ = Sn, An are respectively the totally symmetric 
and the totally antisymmetric representations of U{Ni + A'^2) with n boxes. It will turn out that 
the relevant limit in this theory is slightly different from (8.20) and is given by fixing 

v = r]- = ^^^^. (8.23 

' /c 27ri ^ ^ 
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Positive V will correspond to symmetric representations and negative v to antisymmetric ones. 
In the 't Hooft limit, for fixed N jk^ the two scalings are clearly equivalent. 

The calculation of (8.22) is very similar to the calculation of partition functions of n bosons 
or fermions in the canonical ensemble, where n is fixed and large. But at large n, in the ther- 
modynamic limit, this calculation can be done as well in the grand canonical ensemble. We then 
introduce the fugacity z and consider the grand-canonical partition function, using the expression 
for the determinant as the generating function of the characters 



j; z-W:! = ( det (1 - r^z Ur> ) = ( exp ( J] '^v'^'z' 



n>0 



^.24) 



£>1 



The average value of n in this ensemble is given by (we remove the average notation here, as is 
standard in the grand canonical formalism) 



oz 
This is inverted to determine the fugacity as a function of the number of particles 

^^ — -2-^ yll J ^ 

and then the original VEV can be calculated, in a saddle point approximation, as 

Trf/^ 



Wl « z-"H,(z,) 



exp I - n log z* + 2J 



-r'zi 



^.25) 



^.26) 



^.27) 



For convenience, let us henceforth absorb Y = r/z. It can be seen that, at leading order in 
large iV, the grand-canonical partition function (8.24) is given by disconnected planar graphs. 
Therefore 

(TV V% . 



H,(y)«exp(^5(n), 9{Y)=9sY. 



i 



y^ 



(8.28) 



where the subscript refers to the planar part. We now observe that the function g{Y) is related 
to the planar resolvent in the lens space matrix model (2.8) and (2.10) by 



Y^9iY)=l{MY-')-t) 



log 



1 r 



^(y + 6)(y + l/b) + y/(Y - a){Y - 1/a) 



(8.29) 



Note that compared to uq in (2.10), the sign between the two square roots is reversed. Integrating 
this equation we get 



9iY) 



— log ( - [y^{Y' + b){Y' + 1/6) + V(y' - a){Y' - 1/a) 



(8.30) 



The initial point of integration is chosen to be y = 0, since around that point the integrand 
approaches a constant (^/2 + 0(Y). This guaranties that for small Y the result of the integration 
will be proportional to the 1/2 BPS Wilson loop (2.35). 
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The saddle point equation (8.25) determining the mean value of n is then given by 

i.e., (8.29) 

1 



-2ttw 

2 L 



y/{Y, + b){Y, + 1/b) + V(^* - a)(y* - l/a)] , Y, = rjz,. (8.32) 



This can be solved explicitly in terms of /3, (" or alternatively in terms of B and k. The solution 
reads 

_ iK e-'''(2''+^) / r 16 sin(27ri^) sin(27r(i^ + B))\ 



4sin(27r(i/ + B)) r V^ k"^ J' ^^'^^^ 

The choice of sign is such that y* = when u = 0. We will write 

W;i « exp (A^/5.) (8.34) 

where A^i, which is identified with the action of a brane probe in the large N string/M-theory 
dual, is given by 

A^ = -27rir/z/log(??y) + r]g{Y,). (8.35) 

In the original variables, in terms of ojq, the integral (8.30) is from infinity to a finite position 
Y^^, and represents the effect of adding a single eigenvalue to the system. This fits with the 
standard dictionary [71] identifying a brane with a single eigenvalue. 

This integral gives an expression for the action of the giant Wilson loop, in the limit (8.20) 
which is exact as a function of the 't Hooft couplings. The derivatives of this integral with respect 
to /3 and C can be evaluated in closed form, as in (2.16), in terms of incomplete elliptic integrals. 
The resulting expression can then be studied at the different limits of the ABJM theory as done 
for other observables in earlier sections. 

If we go to the conifold limit, setting A2 = 0, we get an expression for the giant Wilson loop 
in Chern-Simons theory on S^. In that case there exists an exact expression for the Wilson loop 
for all n. As we show in Appendix B, the above derivation in this limit indeed reproduces the 
CS answer. 

We will now discuss the expansion of the result for the giant Wilson loop for large k, since 
this is the strong coupling limit in which one makes contact with the AdS geometry [33]. In 
terms of B and k, the integral (8.30) reads 



f* dY' 



"^^ dY' /I , 

1 + Y'y -e'^'^Y'iK-iismiTrB)) 



+ \ {I- Y'y - e'^'^Y'iK + 4isin(7r5)) 



(8.36) 



where y* is given in (8.33). 

Expanding y at leading order at large k we get 



y. = 2ie--'(2'^+^)^^^^^^ + OiK-^) = ^ ^ e--^^ + 0(k-2) (8.37) 
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This suggests rescaling Y in the integral (8.36) by k, which allows for a systematic expansion in 
powers of K~^. At leading order the integral becomes 

giY,) = - y * ^ (log Vl - e-'^KY' + 0(k-1)) . (8.38) 

This yields 

g{Y,) = ^Li2(e^'^Ky,) + 0{k-^) = ^Li2(l - e""^''^) + 0{k'^) (8.39) 

Another way to get this estimate is to notice that the highest powers of (" in the series expansion 
in y of g{y) are captured by 

g{y) = hi2{Cy) + --- . (8.40) 

Using the dilogarithem identity (B.3) we conclude that the action (8.35), written in terms of the 
original variable n, is 

^-A, = n7rV2X+'^i2B - 1 + ^) + ^ f^l! - Lisfe-^-"/'^)^ + ©(A-^/^.e-^-^^). (8.41) 



Qs ' 2 ' " 47ri V 6 

Notice that this formula does not display the exchange symmetry n ^ N — n for the anti- 
symmetric case r\ = —\. This is because this symmetry is not present for the antisymmetric 
super- representation, as pointed out in [72]. 

The leading order in A in (8.41) is as expected, i.e., n times the action of the fundamental 
string (and n times an extra framing factor). The non-trivial dependence on v only appears at 
subleading order in A, and therefore will not be visible in the supergravity approximation. As 
mentioned above, there are no known 1/2 BPS brane solutions carrying less than /c/2 units of 
electric charge other than fundamental strings. So we expect that the above action describes the 
interaction of these coincident strings. 

For n a multiple of fc/2 (or of fe, if it is odd), we see from (8.33) that Kk = and the integral 
(8.30) is over a full cycle. The argument of the dilogarithm in (8.41) is unity, canceling the 7r^/6 
term. Since Y* passed through one of the cuts C\ or C2, it is now on a different sheet, and exactly 
at the branch point of the logarithm in liJq{^~^\ This happens exactly for the value of n where 
the strings describing the Wilson loop can be replaced by D2-branes, which are the string theory 
incarnation of the vortex loop operators [18]. This suggests that the vortex loop operators are 
related to eigenvalues along the logarithmic branch-cut. It is possible to use our formalism to 
calculate the perturbative and instanton corrections to the these configurations and it would be 
interesting to understand further their significance in the matrix model. 
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A. Normalization of the ABJM matrix model 

Here we shall fix the overall normalization of the matrix model. As explained in the beginning 
of Section 2, to fix the normalization we must fix the coefficient of the cosh in the denominator. 
This term appears as a consequence of integrating out the matter hypermultiplets at one-loop. 
For general supersymmetric Chern-Simons-matter theories, the contribution of a hypermultiplet 
in representation R is given by [4] 

p n=l ^ ' I \ ^ / 

where p are the weights of the representation, and a is the element in the Cartan algebra given 
by 

a = — diag(^i,--- ,/XAr,, z^i, • • • ,iyN^) . (A.2) 

In [4] the one-loop determinant is evaluated up to a multiplicative constant, 

Z[a] = Y[ (Ccosh (7rp(a)))~^/2 . (A.3) 

p 

The constant C can be determined by setting a = in (A.l) 

71=1 ^ ' ^ 

This is a divergent constant, but as usual when considering determinants on compact manifolds, 
we can compute it by using ("-function regularization. Let us define 

The regularization of the quantity appearing in (A. 4) is then — C^(0). An elementary calculation 
shows that 

Cz(^) = -(2^-l)C(s) (A.6) 

where ('(■s) is the standard Riemann zeta function. Therefore, 

-CkO) = -!f^ (A^T) 

and C = 2. 

B. Giant Wilson loops in Chern Simons theory 

Chern-Simons theory on S^ is a particular case of the lens space matrix model when 6=1 and 
the second cut collapses to zero size, i.e., ti = t, t2 = 0. It gives the leading behavior of the 
Wilson loop in ABJM theory when A2 ^ Ai, as discussed in Section 6. 
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Here we consider the behavior of the giant Wilson loops, those in high dimensional symmetric 
or antisymmetric representations presented in Section 8.2, in this limit. In this case it is easy to 
calculate explicitly the action (8.35), since the integral 



g(Y) = -J — log(/i(y')), HY) = -[l + Y + ^{l + yy - 4e^Y\ (B.l 



f-Y 
/O 
can be obtained in closed form 

7r2 1 



(B.2) 



9{Y) = y - 2 log'(^(^)) + ^og{h{Y)) (log (l - e-'h{Y)) - log(l - h{Y)) 
- Li2(/i(y)) + U2{e-'h(Y)) - Li2(e-*). 
Here we used the dilogarithm identity 

Li2(l -x) = — - Li2(x) - log(x) log(l - x). (B.3) 

D 

The solution of the saddle point equation (8.25) is obtained by setting in (8.33) 

K = -4isinh^, B=^^ + l (B.4) 

and we find 

1 — e~'^^^'^ 

The action (8.35) is 

r]Ar, = -2mu log{r]Y^) + g{Y^) 

^2 (B.6) 

= -2mulogr] - 2^'^v'^ + 27rii/t + y + Li2 (e^'''''"*) - l^i2{(?^''') - Li2(e-*). 

Notice that this expression is exact in t. 

We can test (B.6) in all details against a direct calculation of correlators. Indeed, the VEVs 
(Tr/j U) for the Chern-Simons matrix model on S^ are proportional to quantum dimensions (see 
for example [21]): 

{Ttr U) = g''«/2+^(«)^/2dim,(i?). (B.7) 

In this equation, 

q = e^% (B.8) 

ll{R) is the number of boxes in i?, and k/j is the framing factor, given by 

KR = Y,k{k-2i + l), (B.9) 

i 

where li are the lenghts of the rows in the diagrams. The quantum dimensions of the symmetric 
and antisymmetric representations are given by 

»7n(n-l)/4 nt/2 " 

d\m,{ni) = ^ —^ n(l - e-*^-''^^-'^ (B.IO) 



[n]! 



i=l 
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where 

n n 

[n]\ = J](gV2 _ ^-V2) = q\n{n+l) "Q^^ _ ^-^y (g_^l) 

4=1 i=l 

At large n we rescale 

C = -, g-* = exp(-5sO ^ e-^-i'?-? (B.12) 

n 

so that 

log([n]!) « - (-Ti'^v'^ + 2TTir]v f deiog(l - g-^^''''^^)') . (B.13) 

This gives the foUowing contribution to the action 

7r2jy2 + !L _ Li2(e-2'^''''^) = f] (tt^u^ - 2TTw\ogr] + — - Lhie-^''''')^ . (B.14) 

To derive the expression on the right hand side we used, for ry = — 1 the difogarithm identity 



vr^ x^ 



Li2(e") = -Li2(e'") + — - — ± mx. (B.15) 

The product in the numerator of both the symmetric and antisymmetric representations can 
be written in a unified form as 

27rir]iy [ deiog(l - e'^e-'^'''"'^) = r, {Lhie-'-^'"'"') - Li2(e-*)) . (B.16) 

Jo 

The prefactors in (B.7) and (B.IO) contribute 

r/(-37rV2 + 27rizyt). (B.17) 

Together with (B.14) and (B.16) this exactly reproduces (B.6). 

In the antisymmetric representation the result can also be written as 

-27riz/(t + 27riz/) + ^ + Li2(e-*) - Li2(e-*-2--) - Li2(e2--). (B.18) 



This expression agrees at leading order with the D6-brane calculation (8.21) and should be the 
full answer in the limit of A2 = 0. In this expression we see the expected symmetry [68] 

n^N-n (B.19) 

which is 



27riz^ o -t - 27rii/. (B.20) 
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